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Abstract 

We study the constraints of crossing symmetry and unitarity for conformal field theories in 
the presence of a boundary, with a focus on the Ising model in various dimensions. We show 
that an analytic approach to the bootstrap is feasible for free-field theory and at one loop 
in the epsilon expansion, but more generally one has to resort to numerical methods. Using 
the recently developed linear programming techniques we find several interesting bounds for 
operator dimensions and OPE coefficients and comment on their physical relevance. We also 
show that the "boundary bootstrap" can be easily applied to correlation functions of tensorial 
operators and study the stress tensor as an example. In the appendices we present conformal 
block decompositions of a variety of physically interesting correlation functions. 
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1 Introduction 

The "bootstrap" has been a recurring dream in theoretical physics. It is the ambitious aspiration 
that, starting from a few basic spectral assumptions, symmetries and general consistency require- 
ments (such as unitarity and crossing) will be powerful enough to fix the form of the theory, with no 
reference to a Lagrangian. The dual models of the strong interactions emerged as an incarnation of 
the S-matrix bootstrap attempts of the 1960s and eventually led to the discovery of string theory. 
The bootstrap program for conformal field theories (CFTs) in d dimensions was formulated in the 
early 1970s pp. Despite important formal developments such as the operator product expansion 
and the conformal block decomposition (see e.g. the early books [21 [3]), attempts to solve CFTs 
in arbitrary dimensions were not successful. For two-dimensional CFTs, the revolution came in 
the 1980s with the discovery of many exactly-solvable "rational" models. While this is a beautiful 
incarnation of the bootstrap idea, the methods that work in 2d rational CFTdj are too specialized 
to be imitated in higher dimensions, or even in two dimensions for the generic non-rational model. 

The interest in CFT in various dimensions is nowadays stronger than ever, sustained by phe- 
nomenological questions in condensed matter physics (d = 3) and particle physics (d = 4), as well 
as by more formal motivations such as the AdS/CFT correspondence and the rich integrability 
structures of superconformal field theories (d < 6). A pioneering work [4] has rekindled the confor- 
mal bootstrap, turning it into a concrete computational tool. This approach has been refined and 
extended in a series of papers [3ElEllBllSlliniinilI21II3]. 

The modern bootstrap starts with the simple question: in a generic theory, which values of 
operator dimensions and OPE coefficients are compatible with the constraints of crossing symmetry 
and unitarity for the four-point functions? There is a shift of viewpoint, from trying to find analytic 



4 or in closely-related models such as Liouville theory 



2 



answers in a specific model to deriving (by numerical methods if necessary) universal bounds valid 
for any model. As it turns out, one can derive strong constraints already from the analysis of a 
single four-point function of identical scalar operators [1] . This should be regarded as the first step 
in a systematic exploration of the space of CFTs. More surprisingly, important theories such as 
the 3d Ising model appear to live at interesting corners of the parameter space, sitting at "kinks" 
of the exclusion curves [5j [121 [13] . So even the solution of some special models in d > 2 may not 
be too far-fetched, after all. 

In its simplest version, the revived conformal bootstrap works as follows. The four-point corre- 
lation function {ty{x\)ip{x2)<p{x$)ip(x4)) of a scalar operator can be written as a sum over conformal 
blocks in two different channels, by taking OPEs in two different limits. The conformal block de- 
compositions in either channel must sum to the same four-point function, giving crossing-symmetry 
relations for the couplings and scaling dimensions. While this was understood long ago, the main 
idea of [I] is that these constraints can be put to good use by taking derivatives of the four-point 
function at symmetric points and applying linear programming techniques to obtain contradictions 
if certain conditions for e.g. the operator spectrum are not met. The prototypical example of a 
constraint that arises in this way is an upper bound for the dimension of the first scalar primary 
ip 2 appearing in the OPE of two tp's. Crossing symmetry and unitarity imply that A„2 < /(A^,) 
for some numerically determined function /(A„). The method admits straightforward extensions 
to bounds on scaling dimensions of tensorial operators, central charges and OPE coefficients If] 

In this paper we extend this program to conformal field theories with a boundary. An Euclidean 
CFT in d dimensions can be defined in the half-space x^ > 0, with boundary conditions at Xd = 
that preserve an SO(d, 1) subgroup of the original SO(d + 1,1) conformal symmetry [191 I2D]- 
For a given bulk CFT, different consistent boundary conditions are usually possible. Boundary 
CFTs (BCFTs) are very interesting in their own right and find diverse physical applications. They 
describe surface phenomena in systems near criticality, with surface critical exponents related to 
the conformal dimensions of the boundary operators. In string theory, two-dimensional worldsheet 
BCFTs are interpreted as D-branes. These would be sufficient reasons to consider the boundary 
bootstrap, but one of the main questions we would like to address is whether by probing the theory 
with a boundary one can constrain the original bulk theory itselfQ One could in fact also go ahead 

5 Analogous "sum rule" techniques can also be used to obtain non-trivial bounds from modular invariant partition 

functions, see [T4l IPol [TfJl [Tfl fT8] . 

6 A prototype is the the beautiful theory developed by Cardy [211 I22| in 2d rational CFTs, which relates the set of 

consistent boundary conditions with the bulk spectrum and its modular transformation properties. 
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and consider a more general setup where conformal defects of all possible codimensions (boundaries 
being the special case of codimension one) appear on a democratic footing. 

Besides the spectrum of bulk operators and their three-point functions, which are unaffected 
by the boundary conditions, a BCFT is characterized by additional boundary data: the spectrum 
of boundary operators, their three-point functions, and the bulk-boundary two-point functions. A 
correlator containing both bulk and boundary operators can be decomposed in different channels, 
giving crossing-symmetry constraints that in general involve both bulk and boundary data. We will 
focus on the simplest non-trivial type of correlator, the two-point function of two bulk operators, 
which in the presence of a boundary is a non-trivial function of a single conformal cross-ratio. It 
can be decomposed in the bulk channel, by first fusing the two bulk operators together, or in the 
boundary channel, by taking the boundary OPE of each bulk operator. See figured] on page[H 

Outline 

The main advantage of using the boundary bootstrap to constrain bulk dynamics is the simplicity 
of the setup just described. This follows from the results of section [21 where we discuss the two- 
point function of bulk scalar operators: its functional form and its conformal block decomposition 
in the bulk and boundary channels. The conformal blocks turn out to be simple (hypergeometric) 
functions of the single cross-ratio and furthermore depend analytically on the spacetime dimension 
d. This is to be contrasted with the standard conformal blocks for four-point functions (in a theory 
with no boundary) , which depend on two cross-ratios and admit closed- form expressions only when 
d is an even integer. 

In section [3] we demonstrate a remarkable simplification of the boundary bootstrap in a few 
special cases, where one can explicitly solve the bootstrap equations by making an ansatz containing 
only a few conformal blocks in either channel. By this route we are able to recover one-loop results in 
the epsilon expansion purely from the bootstrap equations. These methods do not straightforwardly 
extend to higher loops but give a nice pedagogical illustration of the constraining power of crossing 
symmetry. 

In section |4] we apply the linear programming techniques of [4j to the boundary crossing sym- 
metry equations for scalar two-point functions. We derive a number of general bounds on operator 
dimensions and OPE coefficients. Our bounds however come with a major caveat: while unitarity 
guarantees that the coefficients in the boundary conformal block expansion are positive (since they 
are squares of real numbers, as in [1]), this is not automatically the case for the bulk expansion. 
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Indeed it is not difficult to find counterexamples for certain choices of boundary conditions. We 
need then to assume the existence of boundary conditions where the coefficients multiplying the 
bulk conformal blocks are positive. We present circumstantial evidence for this assumption in the 
appendices where we show that it holds in a large number of calculable cases, for favorable choices 
of the boundary conditions (the so-called "extraordinary" and "special" transitions). It would 
however be more satisfactory to find a general proof. 

External tensorial bulk operators are also more easily incorporated in the boundary setup. We 
illustrate this in the second half of the paper, where we consider the two-point function of two bulk 
stress tensors. In section [5] we discuss the different tensorial structures, the bulk and boundary 
conformal block decompositions and the resulting crossing symmetry equations. We then apply the 
linear programming techniques in section [6] and derive interesting bounds. As before, these results 
rely on certain positivity assumptions for the coefficients of the bulk conformal blocks. 

In appendix O we present a brief derivation of the conformal blocks for a scalar two-point 
function. The remaining two appendices are dedicated to a discussion of a large number of solutions 
to the crossing symmetry equations: we consider scalar two-point functions in appendix [B] and 
stress-tensor two-point functions in appendix O These solutions offer partial justification of our 
positivity assumptions in sections H] and El We also consider an interesting two-point function in 
Liouville theory (with ZZ boundary conditions) that interpolates between all the minimal models. 
We discuss how the analogous bulk four-point function helps to explain a few features of the "kinks" 
observed in the bulk results of [121 02] • 

2 Boundary crossing symmetry for scalars 

In this section we introduce the general setup of boundary CFT and derive the crossing symmetry 
equations for the two-point function of bulk scalar operators. For background material on BCFTs 
see [19j [22], E21 [2H EH], and especially the paper by McAvity and Osborn [26], whose results we 
borrow at several points in this and subsequent sections. 

2.1 Scalar two-point function 

Let us start by deriving the form of the scalar two-point function in the presence of a boundary, a 
classic result dating back to [19]. We will use standard Euclidean coordinates x^ 1 = (x 1 , . . . , x d ) and 
consider the half-space defined by x d > 0, the coordinates tangential to the boundary are denoted 
x. It will be useful to embed this physical space in a higher dimensional space as the so-called 
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null projective cone |27l I28j . Consider Minkowski space in d+2 dimensions in lightcone coordiates 
denoted by P A = (P + , P~ , P 1 , ■ ■ ■ P d ). The null projective cone is denned as, 



P A P A = with P A ~ \P A . (2.1) 
The map from the null projective cone to our physical space is given by 

= p+ " (2-2) 

One easily finds that the usual SO(d + 1, 1) Lorentz group of the d+2-dimensional Minkowski 
space becomes the conformal group of the (i-dimensional Euclidean space. The null projective cone 
provides a linearization of the action of the conformal group. 

As we mentioned above, the presence of a boundary at x d = breaks the symmetry group to 
SO(d, 1). In the null projective cone this breaking can be implemented by introducting a fixed 
vector V with components 

V A = (0,...,0,1), (2.3) 

and restricting ourselves to those Lorentz transformations that leave V A invariant. The residual 
conformal transformations for the coordinates a;' 4 are easily obtained from the linear transformations 
of the P A coordinates. 

Let us now consider scalar fields that are homogeneous functions of the coordinates, 

0{XP) = \- a O(P) , (2.4) 

where A is the conformal dimension of the field O. The physical CFT scalar operator is defined as 

0{x) = (P + ) a O(P) . (2.5) 

The two-point function of O should be invariant under SO(d, 1) and consistent with (|2.4p . The 
only SO(d, 1) invariants that can be formed with two coordinates and the fixed vector V A are 

Pi-P 2 , V-Pl, and VP 2 . (2.6) 

The two-point function must then be of the form 

(Om)(hW) = {2V . Pi) a! {2V . P2) aJ (0, (2.7) 

where /(£) is an arbitrary function of the conformal invariant, 

_ -Pi ■ P 2 

? 2{y ■ Px){v ■ p 2 y K ' } 

6 



In physical coordinates, 

We see that the limit £ — >■ corresponds to bringing the operators close together while the limit 
£ — > oo amounts to bringing the operators close to the boundary. It will be useful to introduce a 
function = £^ 1+ ^ 2 ^ 2 /(£), the two-point function then becomes 

«M«oo.<«.» = p^r'^ff). (2.io) 

For two identical (canonically normalized) operators lining G(£) = 1, since we need to recover the 
usual two-point function far away from the boundary. Although using the null projective cone is 
somewhat of an overkill for the scalar two-point function, this formalism will become essential for 
the tensor calculations of section 

2.2 The boundary bootstrap 

Much like a four-point function for a CFT without a boundary, one can decompose the correlation 
function (|2.10p into conformal blocks. In this case there exist two different decompositions (or 
channels) and we review both of them below. 

In the bulk channel we simply substitute the bulk OPE in the two-point function (|2.10p . For 
two identical scalar operators the bulk OPE takes the form (omitting tensor indices for simplicity): 

0(x)0(y) = 1 + A fc C[x - y, d y ]O k {y) , (2.11) 
{x y) k 

where k labels conformal primary fields. The differential operators C[x — y,d y ] are determined by 
the (bulk) conformal symmetry and the couplings can be taken to be real [3]. We emphasize 
that this OPE is a local property of the bulk CFT and therefore unaffected by the presence of a 
boundary. On the other hand, whereas in the absence of any boundaries only the identity operator 
gets a non-zero one-point function (and all other terms in the OPE therefore drop out of the two- 
point function of O), this is no longer the case once a boundary is present. Using the null projective 
cone it is easily demonstrated that boundary conformal invariance allows for one-point functions 
of scalar operators of the form: 

{ ° {x)) = J0p ' (2 - 12) 

with a coefficient ao whose magnitude is unambiguous as we have normalized the operator using 
the first term in (|2.1ip . One-point functions for operators with spin are not allowed by conformal 
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invariance, see section 15.21 below. Substituting now (|2.1ip in (|2.10p and using (|2.12p one arrives at 
the bulk channel conformal block decomposition: 

G (0 = l + J> fe a fe / bulk (A fc ;£), (2-13) 

k 

where the bulk conformal blocks /bulk(Afc>£) can be determined by working out the expression: 

C[x — y, d y ]j^y^ ■ (2.14) 

This computation was performed in [26], with the result that (see appendix [X] for a new derivation) 

/b u i k (A fc ;0=e Afc/2 2Fi(^,^;A fc + l-^;-e) . (2.15) 

Equations (|2.13p with the explicit expression (|2,15p summarize the bulk block decomposition of 
the two-point function. Notice that the blocks are naturally defined as a series expansion around 
£ = 0, which is when the two operators approach each other. Convergence of the OPE away 
from the boundary however implies that the conformal block decomposition should converge for all 
physical values of £, that is for all < £ < oo. 

In the boundary channel we use the bulk-to-boundary OPE where a bulk operator is written as 
an infinite sum over boundary operators. For a scalar operator this OPE takes the form: 

°( x ) = T^Ka +^iD[a^,dg]6 l (g) , (2.16) 
(zx ) i 

where the index I runs over boundary primary fields, the differential operators D[x d ,dg] are again 
completely determined by (boundary) conformal symmetry and the couplings m are again assumed 
to be real. The first term in (|2.16p corresponds to the one-point function of 0(x) and represents 
the contribution of the boundary identity operator. Subsequent operators all have to be scalars by 
boundary Lorentz invariance. Notice also that in equation (I2.16P we used a hat to denote operators 
living on the boundary (and such operators obviously can depend only on x). 

The constraints of boundary conformal invariance for the correlation functions of boundary 
operators 0(x) are exactly the same as those of ordinary conformal invariance in d — 1 dimensions. 
This implies in particular that boundary operators cannot get one-point functions and their two- 
point functions take the canonical form, 



which also provides a normalization for boundary operators. Combining now (|2.16j) and (|2.10p and 
using (|2.17p one arrives at the boundary channel conformal block decomposition: 



(2.18) 



where the boundary conformal blocks /bd y (A;;£) can now be determined from: 



D[x d ,dg]D[y d ,dy 



1 



(2.19) 



y \x — y| 2A 

Just as for the bulk blocks, this computation was done in [26] (and rederived in appendix IA"1). 

/ bdy (A;e)=r A 2i ? i(A,A + l-^;2A + 2-d;-^ . (2.20) 

The boundary blocks have a good series expansion when both operators approach the boundary, 
that is around £ = oo. 

The boundary block decomposition is summarized by equations ()2. 18j) and (|2.20p . The conver- 
gence of the bulk-boundary OPE away from other operator insertions implies that this conformal 
block decomposition should converge for all < £ < oo as well. 

The statement of crossing symmetry is nothing more than the fact that the two decompositions 
(gJH) and (jXTg]) should agree, 



(2.21) 



A pictorial representation of this equation is shown in figure [TJ The aim of this paper is to explore 
how equation (I2.21|) can be used to constrain the space of boundary conformal field theories. 




E 



EI 1 



Figure 1: Two-point function crossing symmetry in boundary CFT. 
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3 The boundary bootstrap in the epsilon expansion 



In this section we demonstrate that in a few special cases it is possible to obtain an analytic 
solution of the crossing symmetry equation (|2.2ip . As we will see below, in this way we can in 
fact bootstrap the outcome of a one-loop computation and recover the order e critical exponents 
of the Wilson-Fisher fixed point! This is possible because our solutions turn out to have only one 
or two blocks in either channel and equation (12.2ip reduces to a finite-dimensional linear system. 
This should be constrasted with the conformal block decomposition for the bulk four-point function, 
whose asymptotic properties dictate that it always decomposes into an infinite number of conformal 
blocks [4J, which makes the problem much harder. The results in this section therefore highlight the 
relative simplicity of the boundary bootstrap program. At higher orders in the epsilon expansion, 
the problem becomes infinite-dimensional even in the boundary case, and more powerful methods 
will have to be developed. 

3.1 The simplest bootstrap 

Let us begin our exploration of the constraining power of the crossing symmetry equation (|2.2ip 
by considering the following question: is it possible to satisfy crossing symmetry with just a single 
block in either channel? It turns out that this question can be answered affirmatively and leads 
to a rederivation of the free-field theory two-point functions. In formulas, our question becomes 
whether there exists a solution to the equation 

1 + Xa v f hulk {m = C A {ah + fJ 2 /bdyM 0) , (3-1) 

for all £ and with unknowns Aa^, r], A, o? and rj . We use r\ and rj to denote the dimensions of the 
single bulk and boundary operator, respectively. 

In order to find a solution we will below expand both sides in £. The bulk conformal blocks 
(|2.15j) have a natural series expansion in powers of £ around £ = 0, which is when we bring the two 
points close together. On the other hand, the boundary conformal blocks of equation (|2.20p are 
naturally defined via a series expansion around £ = oo where both points approach the boundary. 

Now, using standard hyper geometric transformation formulas (see for example [29J), we can 
expand a boundary block around £ = 0, 

/b dy (V; = ci(i + • • •) + c 2 e~ d,2 H + ...), (3.2) 

with the dots representing subleading integer powers of £ and c\ and C2 certain constants. Substi- 
tuting this expansion into (|3.ip and simply matching the powers of £ to those possibly appearing 
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on the left hand side of (|3.ip . we directly find that: 

A = A^=--1, r J = 2A (/> = d-2. (3.3) 

This is our first non-trivial result: the scaling dimension A has to be that of a free field <f> and the 
value of r\ reflects the simple free-field bulk OPE, x = 1 + </> 2 . 

Our next step is to notice that the bulk block with ij = 2A^ becomes particularly simple, 

/ bu l k (2A^;0= (j^j) \ (3.4) 



and expanding now both sides of (|3.ip around £ = oo we find that 

1 + Aa, + + •••)= ^ (4 + mV (l - I + -..)) , (3.5) 

which allows us to solve for all the other coefficients. We find two possible solutions: 
+ : \a v = +l, Oq = 0, t/ = A0, /i 2 = 2, 
-: Xa v = -1, a 2 o = 0, rf = A^ + 1 , fi 2 = ^— ^ . 

Although we have only used the series expansions of the conformal blocks around the endpoints 
£ = and £ = 00, it turns out that for the above values of the coefficients the crossing symmetry 
equation is miraculously satisfied at every order in £. Therefore, the two functions 

A,,, 



(3.6) 



G + (0 = 1 + /buik(2A^;£) = £ A *(2/ bdy (A ;£)) = 1 + (j^j 
G- (£) = 1- /buik(2A ; £) = ^ Ad y ( A^ + !;£))=!- (-^-) 



(3.7) 



are valid solutions to the crossing symmetry equation (I2.2ip with just a single block in each channel. 
Using (12. lOH we find that they correspond to two-point functions of the form: 

{ftxMy)) = ± , (3.8) 

[x — y) * [x — y r y^<t> 

where y r is the coordinate vector y reflected in the boundary, so if y = (y,y d ) then y r = (y, —y d ). 
This equation informs us that we have derived the two possible two-point functions of a free field 
on a half-space, with the + sign corresponding to Neumann boundary conditions and the — sign 
corresponding to Dirichlet boundary conditions. 

Let us offer a few more comments on the above solutions. First of all, the bulk-to-boundary 
OPE is consistent with the boundary conditions. Indeed, the bulk-to-boundary OPE of a free field 
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4> contains a priori a boundary field <f> and its normal derivative dd<ft of dimensions A^ and + 1, 
respectively. (Notice that these are both SO(d, 1) primaries.) As expected, in the Dirichlet case 
the operator <f> vanishes by the boundary conditions and only the block corresponding to dd<j) is 
present. In the Neumann case the situation is reversed. Finally, the operator 2 is the only operator 
appearing in the bulk channel and the sign of its one-point function is reversed between the two 
boundary conditions. 

3.2 Order e bootstrap 

Having obtained the scalar two-point function for the free theory, let us apply the bootstrap tech- 
nique to the interacting theory in the epsilon expansion. In this section we will allow for N massless 
scalars with strength ^{4> 2 ) 2 ■ The iV-dependence of the free two-point function comes from the 
overall normalization, so the results of the previous section remain unchanged. Defining d = 4 — e, 
the Wilson-Fisher fixed point is given by 

A* 36 „ , 2 



+ 0{e z ). (3.9) 



16tt 2 N + 8 

We can now write the bootstrap equations as a perturbation series in e. Following the strategy 
used in the free case we will assume a finite number of blocks in each channel. In particular, we 
will consider two non-trivial blocks in the bulk channel and a single block in the boundary channel. 
This ansatz has some partial justification in Feynman diagrams. In order for an operator O to 
appear in the bulk OPE of <j> with itself, the three-point function (4>4>0) should be non-zero. For 
operators of the form cj) 2n (ignoring O(N) indices) the only allowed possibilities at order e are (j) 2 
and (j) 4 . For n > 2 the correlator is higher order in e, two or more vertices are needed to contract 
all the legs. In the boundary channe we are only considering the operator 0, similarly to the bulk 
case, the bulk-to-boundary OPE between cf) and </> 2ri+1 for n > is higher order in e. Let us then 
proceed to bootstrap the order e correlator and comment on the validity of our ansatz at the end 
of this section. 

We want to solve the following equation, 

1 + Aa^/ bu i k (A^;£) + Aa^/ bulk (A^; = fx 2 ^f hAy (A^,0 . (3.10) 

Because we are working perturbatively we will write all coefficients as a power series in e. For the 



7 For concreteness we will consider the Neumann case but a parallel analysis can be done for Dirichlet boundary 
conditions. 
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(3.11) 



spacetime dimension d and the external dimension conformal dimension we have 

d = 4-e, 
A = ^-l + SA e + O(e 2 ). 
For the internal conformal dimensions we write, 

A^ 2 = d-2 + <5A 02 e + O(e 2 ), 

A^ 4 =2d-4 + 5A 4e + O(e 2 ), (3.12) 
A^-l + ^ e + 0( e 2 ). 

Finally, for the coefficients multiplying the blocks, 

Aci02 = 1 + SXa^e + 0(e 2 ) , 

Aa 4 =5Xa^e + 0(e 2 ), (3.13) 
M 2 = 2 + 5fi 2 e + 0(e 2 ) , 

where the quantities denoted by "5" correspond to deviations from the free-field solution. For 
example, Aa^4 has only a correction term since it is not present in the free theory. We will again 
use the transformation formulas that led to (|3.2[) in order to expand the boundary blocks around 
£ = 0. The procedure now is the same as before, we Taylor expand both sides of the equation and 
match equal powers of the parameter £. As in the free case, after matching the first few coefficients, 
equation f|3. 10|) is solved to all orders in £. The order e solution is, 

SA <p = 0, 8A.2 = 2a , 5 A, = -a , 

t ( 3 - 14 ) 
5\a<p2 = a , SXa^i = — , 5fi = , 

where a is an arbitrary coefficient. The zero one-loop anomalous dimension for <j) is not a surprise, 
the anomalous dimension of cj) 2 is also well known and can be used to fix the value of a, 

1 /iV + 2\ . 

The first order corrections to the OPE coefficients of the <ft 2 and 4 blocks are positive, while the 
order e correction to fi 2 is zero, as expected from Feynman diagrams. We find a negative anomalous 
dimension for the boundary operator corresponding to (j). The anomalous dimension for (ft 4 does 
not enter the equations at this order in the expansion. The complete corrected two-point function 
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is then 

= * + (jtt) '" + 1 (f+i '^<« + log ^ + >>) + °( £2 > 

= ^(^1(^1) ) A - (2 - £ + £ (£§) ; « + i (£§) W4; « + 0(£2) 

where the + sign indicates Neumann boundary conditions. An analogous calculation can be done 
for the Dirichlet case. We simply quote the result: 

c-u - ' - (^y) 1-1 + k (|±f ) ( " MO + M ? + + <*») 

— (> - ? (&!) ) - « + « (£!) ; « + i (£!) {> + 0(f!) 

= «'~ 4 (( x " 5 + 5 (ft!) ) A*P "5-5 + ° (£2) ' ^ 

which features only minor changes with respect to the previous case. Comparison of these expres- 
sions with the explicit calculation of |26j shows perfect agreement. We have used the bootstrap 
equations to obtain a one-loop result! 

Let us now return to our original ansatz. We did not consider primary operators with derivatives 
acting on the <f>, which we denote schematically by D k (j) 2 and D^ 4 . For the first family, we can never 
have d^d^ acting on the same field, because the equations of motion imply d^d^cp ~ e^ 3 and the op- 
erator is not really of the form d k (f) 2 . The only possibility is to have d^d^ . . . d^cpd^d^ . . . d^cj), 
but these operators are conformal descendants, and their contribution is already taken into account 
by the cj) 2 block. For the second family, the equations of motion argument still holds, but not all op- 
erators are conformal descendants. In fact, there is an infinite number of primaries of the schematic 
form n^ 4 !^ Our original ansatz was thus incomplete, we should have added an infinite number 
of blocks to the left-hand side of equation (|3,10p with tree level dimension = 2{d — 2) + 2k. 
As we obtained the correct answer, it is clear that these operators do not appear at one loop. We 
believe that this is due to the vanishing of the three-point functions (<^O fc ^> 4 ) for k > 0, a fact 
which should follow from the higher-spin Ward identities of the free theory. 

Starting at order e 2 , crossing symmetry can no longer be solved with a finite number of blocks. 
It would be nice to find more powerful analytic techniques to deal with the infinite-dimensional 
linear system, and develop a bootstrap apprach to the all-order epsilon expansion. At each order a 
8 This statement can be checked using conformal characters. 
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new infinite family of bulk primary operators appears. Perhaps the constraints of sligthly broken 
higher-spin symmetry [30|, [3Tj could help in organizing the information contained in (|2.2ip . We 
leave this as an intriguing direction for future work, and devote the rest of the paper to numerical 
investigations. 

Statistical mechanics intermezzo 

In the study of critical systems with a boundary it is well-known that Neumann boundary 
conditions for the Landau-Ginzburg field <f> (which corresponds to the bulk spin operator a) describe 
the so-called special transition, while Dirichlet boundary conditions describe the ordinary transition. 
The phase diagram of the Ising model in the presence of a boundary is shown in figure [2j 
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Figure 2: Phase diagram for the surface critical behavior of the Ising model in dimension 2 < d < 4. 
Temperature is plotted on the horizontal axis and the (relative) surface interaction strength on the 
vertical axis. The extraordinary transition disappears for d = 4, while the special transition is 
absent in d = 2. 

In our investigations the bulk is always critical so we are always on the vertical line in figure [21 In 
the presence of strong boundary interactions the boundary can order at a higher temperature than 
the bulk. The bulk transition where the boundary is already ordered is then called an extraordinary 
transition. For weak boundary interactions one finds the ordinary transition where the boundary 
follows the bulk and also orders at the bulk critical temperature. Finally, there is a critical boundary 
interaction strength where the boundary and bulk critical temperature just coincide which is the 
special transition. There is a third possibility, called the extraordinary transition, which breaks the 
Z2 symmetry of the Ising model, as <j> should acquires a one-point function of the form (|2.12p . The 
extraordinary transition cannot be described in free-field theory (such a one-point function does 
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not satisfy the free equations of motion), but it appears at first order in the Wilson-Fisher fixed 
point in 4 — e dimensions, see appendix IB. 41 We refer the reader to \32\ [33] for introductions to 
boundary critical phenomena. 

The BCFT associated to the extraordinary transition is the most "stable" as there are no 
relevant boundary scalar operators. In fact it is believed that its lowest-dimensional boundary 
scalar is the "displacement operator" T^d, which is the boundary limit of the bulk stress tensor 
with both indices pointing in the direction normal to the boundary. The displacement operator has 
protected conformal dimension exactly equal to d, and it is thus irrelevant on the (d— l)-dimensional 
boundary. The BCFTs associated to the ordinary and special transitions preserve the Z2 symmetry, 
which thus remains a good quantum number for boundary operators. The boundary spectrum of 
the BCFT associated to the ordinary transition contains a single relevant scalar operator which is 
Z2 odd, and corresponds to dd<ft in the Landau-Ginzburg description. Finally there are two relevant 
scalars in the BCFT for the special transition, one Z2 odd and the other Z2 even, corresponding 
respectively to <f> and (j) 2 . 

In d = 2, the extraordinary transition is associated to the Cardy boundary states and |e)) 
labelled by the identity and the energy, respectively. We have 

W), 

(3.18) 

where the kets on the right-hand side denote Ishibashi states. We see that the two states are 
physically equivalent since they are being related by Z2 conjugation. The ordinary transition is 
associated instead to the Cardy boundary state \a)) labelled by the spin, which is given by 

\a)) = \l)-\e). (3.19) 

There is no 2d BCFT associated to the special transition, since the one-dimensional boundary 
cannot order dynamically at non-zero temperature and so the surface transition is absent. 

4 Numerical results for scalars 

Despite the promising results obtained at zeroth and first order in the e expansion, currently no 
good analytic tools are available for the exploration of the general space of solutions of the crossing 
symmetry equation (|2.2ip . Therefore we have to resort to numerical approaches. In this section we 
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adapt the numerical methods of [J] to our case and derive exclusion curves for operator dimensions 
and OPE coefficients. 

The results we obtain below will depend sensitively on some assumptions about the boundary 
operator spectrum and thereby fall naturally into different categories related to the different possible 
boundary conditions. Following [13] we will focus mainly on correlation functions of the a operator 
in the three-dimensional Ising model, whose possible boundary conditions were presented in figure 
[2J For reasons to be discussed in subsection l4.lt our focus will be on the special and extraordinary 
transitions, which will respectively be discussed in subsections 14.21 and 14.31 below. The relevant bulk 
and boundary operator product expansions and scaling dimensions are summarized in table [H For 
d = 4 there are several operators that do not appear in OPE and we indicated this with a dash. 
The quoted values for the Ising model in d = 3 are of course approximate, but good enough for the 
numerical precision of this paper. We were unable to find a reliable estimate of the dimension of 
the a 1 operator for the special transition. 
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Table 1: Bulk and boundary operator product expansions and operator dimensions in the Ising 
model in various dimensions. There is no special transition in two dimensions. For the extraordinary 
transition the first boundary operator is T^d whose dimension is always equal to the spacetime 
dimension d. The results for d = 3 are approximate and were obtained from \13\ [M] whereas the 
results for d = 2 and d = 4 can be found in the appendices. 



4.1 Implementation 

Let us review how to implement the optimization problem numerically. The following techniques 
were explained in great detail in [U [8] so we shall be brief. We start by isolating the contribution 
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of the identity operator in equation (|2.21]h 

1 = " E A ^ /buik(A fc ; e) + C Acxt [a 2 + £ Hi /bdy(A,; £) j , (4.1) 
and introduce the compact notation, 

l = J2p A F A (0, (4.2) 
A 

where 

PA = (\ k a k ,ao,tf) , (4.3) 

F A (0 = (-/buik(A fc ;0^ Acxt ^ Acxt /bd y (A ; ;0) • (4.4) 

With these definitions equation (j4.2j) is analogous to the sum rule of [1] . There is however a crucial 
difference between the boundary problem that we are studying compared to the four-point function 
crossing symmetry of [4] : even assuming unitarity (as we shall always do) the coefficients p A are not 
all guaranteed to be positive. They are certainly positive in the boundary channel, since they are 
squares of real numbers, but in the bulk channel the combination XkO,k is not manifestly positive. 
Indeed it is not difficult to find counterexamples (such as a free scalar with Dirichlet boundary 
conditions). In the following, we will assume positivity for the bulk expansion such that p A > 
as in the four-point function case. The conjecture is that for a given bulk CFT, there exists a 
choice of boundary conditions that exhibits positivity. In the Ising model, the ordinary transition 
is excluded from our analysis, since both signs occur in the bulk expansion (as can be demonstrated 
in d = 2 and in d = 4 — e dimensions). We will however assume positivity for the special and the 
extraordinary transitions. This assumption is supported by the results in the previous section as 
well as in the appendices. We have found positivity of the bulk block coefficients around d = 4, 
both for the free field and the Wilson-Fisher fixed point at order e, as well as in d = 2 where it is 
a consequence of the positivity of the first two coefficients in the first line of (|3. 18j) . In appendix 
IB .71 we also found that the coefficients for the special transition are positive in the O(N) model at 
large N for any dimension. 

We are now ready to start extracting information from the sum rule (|4.2p . The simplest possible 
bound can be obtained as follows: We allow for the bulk spectrum to span all possible values 
consistent with unitarity, 

Abulk > — =- , (4.5) 
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while restricting the boundary spectrum to be greater than a given value, 

A bdy > A min . (4.6) 

Then, we consider a functional A with the following properties, 

A(l) < 0, (4.7) 
A(F A ) > 0, (4.8) 

where, according to our definitions, -Fa stands for any of the blocks appearing in (|4.4p with scaling 
dimensions obeying (j4.5|) and (|4.6p . If such a functional is found, equation (|4.2p becomes incon- 
sistent and we can rule out that particular CFT. The idea then is to see how high we can push 

Before implementing the machinery of linear functionals we need to choose a set of "coordinates" 
in our function space. We will parametrize the blocks by an infinite vector of derivatives {F^} 
evaluated at £ = 1, 



Fa 



(4.9) 



and crossing symmetry becomes now an infinite set of algebraic equations. In order to make the 
problem numerically tractable we will discretize the spectrum of bulk and boundary dimensions 
and consider a maximum number of derivatives. With this truncation we have an optimization 
problem with a finite dimensional set of inequalities, this is an example of a linear program. In 
order to solve the linear programs we used the Mathematica routine LinearProgamming and the 
IBM ILOG CPLEX Optimizer. In all our plots below we used a grid of 5^ = 0.01 and a total of 15 
derivatives. 

4.2 Special transition 

In the following we present our numerical results for the special transition. The one-point function of 
the bulk spin operator a vanishes since the Z2 symmetry is unbroken by the (Neumann) boundary 
conditions. As we have emphasized in the previous subsection, positivity of the bulk channel 
coefficients will be a working assumption. 

4.2.1 Simplest bound for the boundary channel 

Let us start by plotting the simplest possible bound of the form described above. Our only as- 
sumption for the bulk spectrum will be the three-dimensional unitarity bound, Ab u ik > 0.5, but 



19 



otherwise bulk operators of any dimension are allowed to appear in the OPE. Crossing symmetry 
and positivity however imply that the conformal dimension of the lowest dimension boundary op- 
erator cannot be arbitrary. Instead, we found that depending on the external dimension the first 
boundary operator has to lie below the curve of figure El 




0.52 0.54 0.56 0.58 0.60 0.62 0.64 

A ex t 

Figure 3: Upper bound for the first boundary operator in the special transition. 



Although this is a correct bound, we should mention the following caveat: The bulk block blows 
up at the unitarity bound and our more precise assumption for the bulk spectrum was actually 
Abuik > 0.5 + 10~ 6 . Unfortunately, it turns out that the numerics are quite sensitive around this 
point. For example, the bound becomes much stronger if we change our assumptions on the bulk 
spectrum to Abuik > 0.51. Because of this, we do not consider this plot to be physically very 
relevant but it serves as a good warm-up example before tackling the most interesting cases below. 

4.2.2 Improved bound for the boundary channel 

The boundary bound obtained above can be improved by making further assumptions. In the 
bulk channel decomposition of a scalar two-point function we expect, on physical grounds, a "gap" 
between the unitarity bound and the conformal dimension of the first operator appearing in the 
bulk OPE. For example, according to table [H in the three-dimensional Ising model the first bulk 
operator appearing in the OPE of the spin operator a is the energy operator e with A e = 1.41, far 
above the unitarity bound. Clearly, allowing for the bulk spectrum to go all the way down to the 
unitarity bound is very unphysical. In figure S] we present an improved bound in which we assumed 
that the bulk spectrum satisfies Abuik > 2A cxt . 
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Figure 4: Improved bound for the first boundary operator in the special transition. The bulk 
spectrum is assumed to satisfy Abuik > 2A ex t- 



Our solution seems to indicate that the bound cannot go below the straight line where Abdy = 
A ex t. The reason for this is the trivial solution (x\ — X2)~ 2Acxt which we discuss in appendix IB. 51 
This two-point function contains no non-trivial bulk blocks and thus effectively has an infinite gap 
in the bulk spectrum. On the other hand, it also has a boundary channel expansion which starts 
with a block of dimension A ext and our bound of course cannot get past this particular solution. 
In a sense, the bound is optimal in this case, going down until it hits a known solution to crossing 
symmetry. 

For the Ising model the dimension of the first boundary operator has a value of ~ 0.42 and is 
well inside the allowed region of figured Ideally, we would have found a plot with some striking 
feature around this value, like the kink of |13j . However, in our case the trivial solution is standing 
in the way. A qualitative explanation for this difference appears in the epsilon expansion results. 
Namely, the anomalous dimension of the e operator (which is <p 2 in d = 4) is positive at one loop, 
so the Ising model lies above any trivial (mean field-like) solutions for the bulk four-point function. 
On the other hand, the one-loop anomalous dimension of the first boundary operator is negative, so 
we end up below the trivial solution. This was of course largely a coincidence - we are not aware of 
any fundamental reason requiring these anomalous dimensions to have a definite sign. Some effort 
was made in order to circumvent the trivial solution but we did not succeed in obtaining reliable 
"kinks" that highlight the presence of the Ising model. 

We would like to stress however that our plot is still teaching us something very non-trivial: 
the lowest boundary dimension can never be greater than the external dimension. Interestingly, 
this result precisely implies that the bulk-to-boundary OPE is never regular, see equation (I2.2ip . 
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It would be very interesting to find a more direct argument for this result — perhaps even one that 
does not rely on our specific assumptions. 



4.2.3 Bounding the second boundary operator in the Ising model 

Our assumptions in the previous section were almost minimal, and the result is a general bound 
valid on the space of BCFTs. In this section we will take a closer look at the three-dimensional 
Ising model and attempt to bound the second boundary operator. We will do so for both the {era) 
and (ee) correlators. Using the results from tabled] we can assume that 



A ext = 0.518 ; 
Abuik > 1.41 , 

A[ 2 i > A (2) . 

bay — mm 



(4.10) 



In the boundary channel the first block corresponds to a. We assume that it sits isolated at 
Av.j„ ~ 0.42 and that all the subsequent blocks have a scaling dimension greater than A^ . 



"bdy 



(2) 

Proceeding as before we push AV n as high as possible until the CFT becomes inconsistent. This 
will give us an upper bound for the dimension of the second operator a' , only valid for the {era) 
correlator of the 3d Ising model. Because ~ 0.42 is our less precise value we will explore a range 
around this number. Our result is shown in figure [5] 
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Figure 5: Upper bound for the dimension of the second boundary operator in (era) as a function 
of the dimension of the first boundary operator. 



The same can be done for the (ee) correlator. The statistical mechanics data [M] in this case 
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Acxt = 1.41 , 
A bulk > 3.80 , 
AS y ~0.75, 
A[ 2 i > A (2 > . 

bay — mm 



(4.11) 



and the resulting bound is shown in figure 

Unfortunately, we were unable to find reliable estimates of the scaling dimensions of the second 
boundary operators in the statistical mechanics literature. It would of course be interesting to 
compare our values with e.g. a two-loop computation for the Wilson-Fisher fixed point. 
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Figure 6: Upper bound for the second boundary operator in (ee) as a function of the first boundary 
operator. 



4.3 Extraordinary transition 

In the extraordinary transition the boundary identity operator is always present, so bounding the 
lowest boundary dimension is not an interesting exercise in this case. The second boundary scalar 
operator is expected to be T^d, the energy momentum tensor with indices in the normal direction, 
evaluated on the boundary. This operator is always present in the boundary spectrum and has 
conformal dimension exactly equal to d, see [24\ for details. Having so much information about 
the boundary channel we would like to address the following question: can we bound the bulk 
spectrum using the boundary bootstrap? We will show below that this is indeed possible, although 
our bound is weaker than the one obtained in [13j who used the crossing symmetry equations for 
the bulk four-point function. 
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4.3.1 Bound for the bulk channel 



The assumptions for the extraordinary transition are 

Abulk > A m 

(i) 



A 



bdy 

(2) 
bdy 



0. 



(4.12) 



A[ 2 i>d, 



where we used a notation familiar from the previous subsection. The fact that A^ = corresponds 
to the boundary identity operator which sits isolated, and we then allow for any operator with a 
dimension greater than (or equal to) d to be present in the boundary channel. A m ; n is the lowest 
bulk dimension and the quantity we want to bound. In figure [7] we plot our bound as a function of 
the external dimension. 

Abulk 2.4 - 
2.2 - 




Figure 7: Bulk bound for the extraordinary transition as a function of the external dimension. The 
dashed line corresponds to the (stronger) bound obtained in [T3] using the bulk crossing symmetry 
equations. 



Because figure [7J can be directly compared with the bound of [13] we have superimposed their 
result on our plot. We can see that the bound obtained using the boundary bootstrap is qualitatively 
different, it is weaker and has no kink at the Ising point. Since we successfully found an "optimal" 
bound for the boundary spectrum in the previous subsection, it is surprising that our bulk bound 
does not exhibit any of the expected features. 

There are two possible explanations for the discrepancy seen in figure [71 First, there may be a 
spurious solution to crossing symmetry that we have not found yet and that prevents the bound 
from going lower. If such a solution exists then it would be interesting to understand whether it 
corresponds to a full-fledged BCFT or not. Notice that this solution would appear to violate the 
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bound of [13] but this may be due to the fact that certain operators do not get one-point functions 
and therefore do not appear in our bulk block expansion. The second explanation is that our 
numerics are not precise enough and that we would be able to lower the bound by increasing our 
numerical precision. We offer some comments on this second possibility below. 

Bulk bound for arbitrary d 

One of the advantages of studying the boundary problem is that the blocks are an analytic function 
of d. In figure [8] we plot the bulk bound obtained above for different dimensions including non- 
integer values. 




8.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Figure 8: Bulk bound for different spacetime dimensions in the extraordinary transition. We 
highlighted the Ising model in various dimensions with the crosses. The dashed line is a specific 
solution for d = 2 which interpolates through the minimal models, see appendix IB. 21 

The bound we find is always significantly different from any known solutions to crossing symme- 
try. In particular, in the figure we have shown the line interpolating through the minimal models in 
d = 2 and the Ising model for the integral dimensions. Again, it would be interesting to understand 
if this is due to our finite numerical precision or whether there exist 'spurious' solutions to the 
crossing symmetry equations at the current bounds. 

4.3.2 Upper bound for Tdd OPE coefficient 

The method of linear functionals can also be used to bound OPE coefficients. In [6] a universal upper 
bound for the OPE coefficient of three scalars was found using the four-point function bootstrap. 
The same technique was used in [7J [8] to obtain an upper bound for the OPE coefficient of the 
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stress tensor. This coefficient is inversely proportional to the central charge c of the theory so the 
result translates into a lower bound for c. 

In this section we will use the boundary bootstrap to bound the coefficient fj% of the Tdd bound- 
ary block fbdy(d,C)- We recall that this block is always present in the extraordinary transition, see 
the OPE in table Q3 We start by imposing, 

A(C Acxt /bd y (d,e)) = 1, (4.13) 
A(F A ) > 0. (4.14) 

Applying this functional to the crossing symmetry relation (I4.2h we obtain, 

A < A(l) , (4.15) 

where /i^ is the OPE coefficient of fbdy(d, £). The best bound is obtained by minimizing the action 
of A on the identity. For the spectrum we require, 
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(4.16) 



Notice that we have again assumed a gap of 2A cxt in the bulk. We plot our result as a function of 
the external dimension in figure [9l 
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Figure 9: Upper bound for the coefficient of the T^d block as a function of the external dimension. 
The dashed line represents and improved bound with a stronger assumption for the gap, following 
the dashed line of figured] (see text). 



Let us try to justify our choice of Ab u ik > 2A ext . A way to make the bound stronger would 
be to increase the bulk gap above this value, the maximum value we can assume for the gap is 
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dictated by the bulk bound of [13] . obtained using the four-point bootstrap equations. In figure 
[9] we have thus plotted an improved upper bound (dashed line) assuming A bu ik > f(^-ext), where 
/(A) is the function represented by the dashed line of figure It is clear that the upper bound is 
not too sensitive to the assumed gap. For example, for the Ising model A ext = 0.518, and the upper 
bounds are $ < 0.0734 and n d < 0.0693 for A bu i k > 2(0.518) ~ 1.04 and A bulk > /(0.518) = 1.41 
respectively. A change of ~ 0.37 in the bulk gap translates into a change of ~ 0.0041 in the bound, 
so at least for this example 2A ext does a good job as a representative gap for the space of CFTs. 

The procedure used above generalizes with no major changes to arbitrary dimensions, let us then 
make a quick comparison with some known values. For the 2c? Ising model the coefficient /j, d can 
be read from the conformal block expansion in (|B.5p . it has the value fi d = s ^/2 ~ 0-0221 whereas 
the Linear Programming methods result in an upper bound /jL d < 0.0309. For the extraordinary 
transition in the e-expansion equation (1B.39|1 tells us \j? d = = 0.10, whereas we obtained the 
upper bound fjR < 0.119 in four dimensions. We see that the numbers agree reasonably well. 

4.3.3 Towards the Ising model 

In analogy with [13] we may try to isolate the Ising model in various dimensions. To this end we will 
improve the results of the previous subsection by using as additional knowledge the dimension of the 
next scalar operator e 1 which appears in the a x a OPE. According to table [Tj in three dimensions 
this operator has a scaling dimension A e / of approximately 3.84 whereas in two dimensions it 
has dimension 4 (it corresponds to L^L— 2l)- We again assumed a boundary channel spectrum 
consistent with the extraordinary transition, i.e. a possible one-point function and a gap equal to 
the spacetime dimensions d. Summarizing, 

A b l > A, , 

Affi y = 0, (4.17) 

with A e / fixed to the values of table [TJ Our aim is now to find the possible range of values that 
Z^bulk can take. The resulting plots are shown in figure [TUl 

Notice that the plots give results that are qualitatively similar to those of [13] , in a considerably 
simpler setup. This is of course an encouraging result. Furthermore, we also did not rule out the 
Ising model and this provides some a posteriori justification for our assumption of positivity in 
three dimensions. 
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Figure 10: Locating the Ising model in d = 2 (left) and d = 3 (right). The plot show the dimension 
of a bulk operator versus the external dimension. With the assumptions explained in the main 
text, we need at least one bulk operator in the shaded regions. The Ising model is indicated with 
the cross in both plots. 



It is however rather unfortunate that the bounds we obtain are relatively weak. For this specific 
example we have tried different numerical implementations as well, for example we have tried to 
include more derivatives or to evaluate the blocks at different points like £ = 1/2 or £ = 2. In 
each case we were unable to significantly lower the bounds. We have also attempted to improve 
the results by imposing an additional gap between the second and the third operator in the bulk 
channel. The third bulk operator has scaling dimensions 8 in d = 2 and approximately 4.6 in d = 3. 
Imposing this additional gap significantly improved the bounds for d = 2 but unfortunately this 
was not the case for d = 3. 



5 Boundary crossing symmetry for stress tensors 

In section [2] we derived the crossing symmetry equation (|2.21|) for two-point functions of scalar 
operators using the bulk and boundary conformal block decompositions. In this section we will 
derive a similar equation for the two-point function of the stress tensor. We will then use this 
equation in section [6] to obtain numerical bounds for the spectrum of operators appearing in the 
stress tensor OPE. 

The main results of this section are summarized in subsection 15.11 We then present the details 
of our computations in subsections 15.21 through 15.41 These latter subsections are not essential for 
the remainder of the paper and can safely be skipped by the casual reader. 
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5.1 Summary of results 



As we show in equation fl5. lTf) below, the two-point function of a spin two operator in the presence of 
a boundary features three independent tensor structures. Each tensor structure comes multiplied 
with its own scalar function of £ and we find it convenient to collect these three functions in a 
three-component vector of the form (/(£), <?(£), Furthermore, for the stress tensor the Ward 

identities relate the three components in the following way: 



d 



d 



(d - 2)£ 2 — g = (d 2 + 3d - 2)h - 2{d - 1)£(1 + £)—h 



d£ 



d£ 

4(1 + £)h + U{d 2 + 2d - 4) - 2d£ 2 (l + £) 



(5.1) 



g, 



so up to a few integration constants there is effectively only one independent function of £. 

In the following subsections we derive the conformal block decompositions of the functions 
(/, g, h) in the bulk and the boundary channel. The main result of these subsections will be the 
following crossing symmetry equation: 

A\ /7buik(A fe ;£)\ 
(.) +^2^kao k 



M( ) 



5buik(Afc;£) 
\/ibuik(A fc ;£)/ 



//£ v (A n ;£)\ 



(5.2) 



bdy 

<?S y (A n ;£) 

V4 2 i y (A n ;£)7 



where all the functions (/, g, /i) are explicitly known functions of £. Equation (15. 2p is the analogue 
of (|2.2ip for scalars and we will use it in section [6] to obtain bounds on operator dimensions and 
OPE coefficients. Let us now discuss it in a bit more detail. 

First of all, because of the three independent tensor structures we get a three-dimensional vector 
of equations (and the conformal blocks themselves also become three-dimensional vectors). It is 
then important to realize that the Ward identities are operator equations and therefore they must 
be true for the individual conformal blocks as well. Each vector appearing in (15. 2D thus individually 
satisfies the Ward identities (15. ip . 

The left-hand side of (|5,ip is the bulk channel conformal block decomposition. As in (|2.2ip . we 
separated out the conformal block corresponding to the identity operator. For the other operators 
we should recall that SO(d, 1) conformal symmetry dictates that only scalars can get non-zero 
one-point functions and therefore only scalar blocks can contribute to the bulk channel expansion. 
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The right-hand side of (|5,ip represents the boundary channel conformal block decomposition. A 
priori, a spin 2 operator has a boundary OPE decomposition involving operators with spins ranging 
from to 2 and indeed we find all these possibilities in (|5.ip . where the spins of the exchanged 
operator is written as the superscript in parentheses. However in this case the Ward identities turn 
out to further constrain the conformal block decomposition. More specifically, the boundary scalar 
and vector appearing in the boundary OPE decomposition of T^ v must have scaling dimensions 
equal to the spacetime dimensions, so A^ ) = A^ = d. There is thus a unique block for the 
exchange of a scalar of dimension d and also for a vector of dimension d. These two blocks are the 
first two terms on the right-hand side of (15.2p . On the other hand, the dimensions of the spin 2 
fields are not constrained in this way and there can therefore in principle be infinitely many spin 2 
blocks, represented by the final sum in (15. 2p . 

Let us offer a few more comments on the spin and 1 boundary operators. As one may 
have anticipated, in physical theories they correspond to the Tdd and Tid {i being a tangential 
index) components of the bulk stress tensor, restricted to the boundary. These operators are 
intimately related to infinitesimal variations in the location of the boundary surface which explains 
the 'non-renormalization' of their scaling dimensions, see |24j for details. For physical BCFTs the 
displacement operator Tdd is generically present on the boundary and we encountered it already in 
the discussion of the extraordinary transition in section [3l On the other hand, the vector operator is 
only present if there is a non-zero energy flow across the boundary. For BCFTs this is an unphysical 
boundary condition and we can then set fj%^ = 0. (Notice that an energy flow would be allowed 
if the surface x d = was actually an SO(d, 1) preserving interface between two different theories, 
one defined for x d > and the other for x d < 0, and in such cases the vector block will generically 
be present.) 

In appendix Owe present a few explicit solutions to the crossing symmetry equation (|5.2p . We 
discuss the universal solution in two dimensions (which is fully determined by the Virasoro algebra) , 
the free-field theory solutions in d dimensions and the extraordinary transition to leading order for 
the Wilson-Fisher fixed point. 

5.2 Correlation functions of tensor operators 

In this section we discuss correlation functions of operators with spin in conformal field theories. We 
will use the results of |35j . see also |36j . and adapt them to conformal field theories with a boundary. 
Many of the results in this and the next two subsections were also obtained in [26, 24J but we present 
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here an independent derivation which is straightforwardly implemented on a computer. 

The index structures appearing in correlation functions of tensor operators are easily found in 
the null projective cone formalism discussed in section [270 According to [35], a generic tensor field 
fn±...v n { x ) lifts to a tensor field Fa^.^a^P) in the null projective cone with the following properties: 

— equal symmetries in the indices of FA 1 ...A n {P) and of j \i\...v n {%)\ 

— transversality, so P Ai FA 1 ...A i ...A„(P) = for 1 < i < n; 

— a gauge equivalence defined as Fa 1 ...a„(P) ~ PlA^ ^ A for any A and 1 < i < n. 

For symmetric traceless tensors it is convenient to contract the indices on F with auxiliary vari- 
ables Z A and write F(P, Z) = Fa 1 ...a ti {P)Z Ai ■ ■ ■ Z An . Tracelessness implies that we may restrict 
ourselves to the subspace defined by Z 2 = and the gauge equivalence implies that we may take 
Z ■ P = as well. The transversality condition becomes: 

P~F(P,Z) = 0. (5.3) 

Correlation functions of n symmetric traceless tensor primary operators can now be written as 
scalar functions G(Pi, Zi) with 1 < i < n with the following properties: 

— the dependence on Zi should be a homogeneous polynomial of degree U; 

— the dependence on Pi should be homogeneous of degree — Aj; 

— transversality dictates that Pi ■ dzfi = for 1 < i < n; 

— for any conserved tensor there is a Ward identity of the form [35] 

(dp- DW)G = 0, (5.4) 

with 

^-H-^-^)4x-WJ^=. (5.5) 



,2 8Z)dZ A 2 dZ-dZ' 

where P and Z are the variables corresponding to the conserved tensor, for example Pi and 
Z\ if the conserved tensor is the first operator. 

As an example, let us review the well-known result for the three-point function of two stress ten- 
sors and one scalar operator Gtto(Pi> Pi> P31 Zi, ^2)- The first three constraints together dictate 
that there are three different invariant tensor structures, 
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with for now arbitrary constants a, b, c and with building blocks 

Uz x ■ P 2 ){Pi ■ Ps) - (Zi • P 3 )(Pi • P 2 j) {{Z 2 ■ P 1 )(P 2 • P 3 ) - (Z 2 ■ P 3 )(Pj ■ P 2/ 

w 12 v 

H±2 



(Pi-P 2 )(P2-P 3 )(P3-Pi) (5.7) 
(Z 1 -Z 2 )(P 1 -P 2 )-(Z 1 -P 2 )(Z 2 -P 1 ) 



Pi -Pi 

The Ward identities for the stress tensor furthermore dictate that 

_ A(A + 2) 



Xtto , 

Xtto , (5-8) 
Xtto , 



4d(d + l) 

(A - d) 2 (d-l) - 2d 

d{d + l){d-2) 
A((A-d)(d-l)-2) 



d(d + l){d-2) 

where Xtto is an undetermined overall coefficient. Upon sending A-f Owe find that a, c — > but 
b — > Xtti and we recover the unit normalized stress tensor two-point function, 

(T(P 1 ,Z 1 )T(P 2 ,Z 2 )) = { _ 2 p^p 2)d , (5-9) 

provided we set Xtti = 1- The normalization in (|5.8p is therefore such that Xtto is a natural 
three-point coupling coefficient. 

Let us finally take the OPE limit by sending Pi — > P 2 . In that case H\ 2 remains finite whilst 

and therefore 

a^*) 2 + feff 2 2 + cg 12 Wg PE 
^ (_2P 1 . j P 2 )^A/2 ( _ 2 p 1 .p 3) A > I 5 " 11 ) 

and we infer that the TxT->0 operator product expansion becomes to leading order 

T(Pi, Z 1 )T{P 2 , Z 2 ) ~ . . . + — - - p2)d _ A/2 O(Pl) + • • • (5.12) 

where we assumed that O is normalized such that {0{P\)0(P 2 )) = (— 2Pi • -P 2 )~ A - 

As we mentioned in section [27TT the breaking of SO(d + 1, 1) to SO(d, 1) due to the presence of 
a boundary is implemented by introducing an additional fixed vector 

V A = (0,0,..., 0,1), (5.13) 

representing the unit normal to the boundary. Correlation functions are still required to be SO(d + 
1,1) scalars with the same four properties as above but they can now depend on V A as well. For 
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example, we have already mentioned that the one-point function of a scalar operator can take the 
form: 

{ ° {p)) = ly^W ' (5 - 14) 

with arbitrary coefficient ao- For one-point functions of tensor operators one directly sees that 
the numerator would have to involve a factor (V ■ Z) 1 but this is not transverse and so higher-spin 
one-point functions must vanish. 

With two points we can build the invariant object £ of section [27T1 which we recall was 

-Pi ■ g _ (*i ~ ^? , , 

K 2{V ■ P X ){V ■ P 2 ) 4xjx d 2 ' 1 J 

and conformal symmetry thus determines two-point functions only up to arbitrary functions of £. 
For the scalar two-point function this leads to equation (|2.10p which was: 

<0i(Pi)0 2 (P 2 )> = {W . Pi) a\ 2V . P2) aM oM , (5.16) 

where /o 1 o 2 (C) is n °t fixed by conformal symmetry. Two-point functions involving tensors are 
easily found, e.g. 

. _ ({Z 2 -V){P 2 -P l )-{P 2 -V){Z 2 -P l )f 

= ^ {v . Pi) a 0+2{v . P2) a t+2 J for(0 , (g i?) 

fjj(0Hi2 + gjj(0Qi2 



ZfZ 2 B (J A (P 1 )J B (P 2 )) 
ZfzfzCz? (Tab (Pi)T CD (P 2 )) 



^(v -Pi) Ai (y -p 2 ) A2 ' 



(4£)A l( y. Pl) A l( y.p 2 )A 2 
with ffi2 already defined above and with 

Q _ ( o^^_ ( ,. Zi))( (L«_ (F ., 2)) . ( , 18) 

If the above tensors are conserved then we write J and T instead of J and T ■ In that case 
A j = d — 1 and At = d and from the Ward identities we also find that: 

foM) = coM{i + i)r d/ \ 

fo T (o = c mi+or 1 - d/2 , 

= (id + 1) - 2£^) 5JJ - 2£ 2 ^ (/jrj + 5JJ ) , (5 . 19) 
(d - 2)^ 2 g' TT = (d 2 + 3d - 2)h TT - 2(d - 1)£(1 + £)/^ T , 

Mef TT = "4(1 + i)h TT + (t{d 2 + 2d - 4) - 2d£ 2 (l + £)4)<?TT , 
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with c .. denoting an integration constant. We see that the two-point function of two stress tensors 
and the two-point function of two currents are both fixed up to a single function of £. The last 
two equations in (|5.19p were already presented in equation (|5.ip . They agree with equation (2.27) 
and (2.31) of [26j with the replacements /(£) = C(v), g(£) = Av 2 B(v) and h(£) = v 4 A(v) and with 

« 2 = e/(e + i). 

We can also insert operators at boundary points labelled X satisfying X ■ V = 0. As before, 
we will denote such operators with a hat. We project the indices of such operators to lie along the 
boundary, which in the null projective cone is implemented by the constraint V • Z)W =0 with the 
operator DjP already given by (|5.5j) . The correlation functions of interest are those with a single 
stress tensor in the bulk. We find: 

2 



(> (Z 2 -V){P 2 -X 1 )-{P 2 -V){Z 2 -X 1 ) 
Z$Z*{0{X x )T AB {P 2 )) = 5 dA . c 6T ± ( _ 2Xi ; ^ 2 



, L , r . (Z 2 -V)(P 2 -X 1 )-(P 2 -V)(Z 2 -X 1 ))H 12 {52Q) 
Z x Z 2 Z 2 (Ja(Xi)T bc (P2)) = SdAj c jt {-2Xi ■ P 2 ) d+1 ' 

jj2 1_ q2 

z£z?2§Zg(tAB(Xi)TcD{P*)) = cfr { _ 2Xi .p 2) 4\v -P 2 ) d ' A t ' 
with 

^ = (Zx • Z 2 )(P 1 • P 2 ) - (Z, • P 2 ){Z 2 • FQ Z^Zt-{Z 1 -V)V A . (5.21) 

Notice that for scalars and vectors the scaling dimension is required to be d whereas the dimension 
of v is unconstrained by the Ward identity. 

Up to terms that ensure that V -D^ annihilates the correlator, two-point functions of boundary 
operators are of the same form as two-point functions of bulk operators in the absence of a boundary. 
In particular we find that: 

1 



(0(X 1 ,Z 1 )0(X 2 ,Z 2 )) 
(JiX^Z^JiX^Z^) 

(f{X u Z x )f(X 2 ,Z 2 )) 



(-2X! • X 2 )^ ' 

H 12 - {V • Z{){V • Z 2 ) 

(-2Xi • X 2 ) A ' (5.22) 

H l2 - {V ■ Z 1 )(V • Z 2 )) 2 - ^(V ■ ZyfiV ■ Z 2 f 



(-2Jfi • X 2 ) A 

Equation (|5.22p defines our normalization of the boundary operators. Notice that H\ 2 descends 
from the projective cone to z^z 2 (5^ v — 2xi 2t/J ,xi 2> v /x\ 2 ) so it is easily verified that our normalization 
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is consistent with reflection positivity. Using (|5.20p and (|5.22p we find the bulk-to-boundary OPE 
of the stress tensor, 

T(P, Z) -> c 6T (Z ■ V) 2 6(X) - Cj T (Z ■ V)J(X, Z) + {v C ^ Af t(X, Z) + ... (5.23) 
5.3 Bulk channel blocks for the stress tensor 

In this subsection we compute the conformal blocks for the two-point function of the stress tensor 
using the conformal Casimir differential equation method of |37J. These are the conformal blocks 
appearing on the left-hand side of (|5.2|) . 

On a symmetric traceless tensor F(P,Z) the action of an element L AB of SO(d + 1, 1) takes 
the form: 



L AB F(P,Z) = (P AWB - P BWX + -^-—(Z A D%+*> - Z B D^))F(P,Z) , (5.24) 

with the operator D^ +2 ^ given by (|5.5p but with d — >■ d+2 since we are rotating in d+2 dimensions. 
The conformal Casimir equation is then: 

h AB L AB F(P, Z) = -C A>l F(P, Z) , (5.25) 

with Ca,z = A(A — d) + 1(1 + d — 2). We used this equation in appendix [A] to find the result (|A.6P 
for the conformal block in the bulk channel for a scalar two-point function. For two stress tensors 
the conformal block can be written as: 

G b (P 1 ,P 2 ,Z l ,Z 2 )- md(y . Pi)d(y . p 2)d , (5-26) 

and the constraint l(L { j B + L^ B )(L^ AB + L^ AB )G A = -C A ,oG A together with the Ward iden- 
tities leads to the unique solution for the coefficients: 

A(A + 2) 



| (4£) A / 2+2 2 iq ( 2 + 1 , 2 + j; 1 - ~ + A; -Z) , (5.27) 



16d(d + 

with fb and g b determined by the Ward identities (|5.19j) . Let us verify the normalization by taking 
the OPE limit £ — > 0. We already mentioned that H\ 2 then remains finite and it is not hard to 
find that 

Qi2^-^< 2 PE , (5.28) 
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with W$, PE defined in ([57TO]) . Prom the expansion of (|5.27p and the Ward identities we find 

and the entire block behaves as: 

b K p L>ft,Z 1} Z 2 ) - ( _ 2Pi .p 2)d _ A/2(y . Pi)A , (5-30) 

which is compatible with (JESJ, (pTT2j) and (JoTI]) . 

Explicit expressions for /j, and 56 are also available in terms of linear combinations of 2-^1 
hypergeometric functions. 

The identity block can be found by sending A — > 0. We then find that /& = 1 and gy = hf, = 0. 

5.4 Boundary channel blocks for the stress tensor 

We label the boundary block associated to a primary operator of dimension A and spin I as Gg 1 ' 
(with a subscript "s" for surface). Each block has again the same form as the TT two-point function 
given in (]5. 17|) with three associated functions fs^ , gs and hi^' l \ In the two-point function of 
the stress tensor there are three types of boundary blocks, G^' ^, Gs and Gi A ' 2 \ To find these 
blocks we act with the SO(d, 1) Casimir operator on one of the two points and solve the resulting 
differential equation. In the equations below we use h = d/2. 

For a block corresponding to the exchange of a boundary scalar of dimension d we find: 

^' 0) = 2h{2h + + tr h - 3 (2H2h + l)f + 2(2h + l)(h- 1)£ + h(h - 1)) , 

f(d,0) _ 1 ph-lf-i , t\-h-l 

Js Ah(2h + lf { 

where we already fixed the normalization. In the limit where £ — > 00 we find that only the third 
tensor structure contributes and 

G^'°) (p 1 ,p 2 ,zi,z 2 ) ~ { Ym V p 2 r )2 ' (5 - 32) 

which agrees with ()5.23j) and the first equation in (|5.22j) . 
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G^(P 1 ,P 2 ,Z 1 ,Z 2 ) v nn nN2fe J -, (5.34) 



For the block corresponding to the exchange of a boundary vector of dimension d we find: 
^ = 2(2fr 1 +l) ^ +1(1 + ^ ^ ( " 2(2/1 + ^ + 2h[k ~ 1)e + Kh ~ l) 

9 ^ 1] = (2/TTi) + °' k ' 2 + h{l + 2e + 2?2) ) ' (5 - 33) 

= 4 (2 ^ +1) ^" 1 ( 1 + 0"* _1 (1 + 2 - 
and the block behaves for £ — > oo as 

(V ■ Zi)(V • Z 2 )[Hx2 - (V • Z 1 )(V • Z 2 ) 
(-2Pi • PV) 

which is again consistent with the formulas given above. 
Finally, for the spin two blocks: 

^ A)2) = 2(h - 1) ( 4 ^ ) 2ft-A 3F2 ^2 + A, 3 - 2/t + A, 1 - /t + A; 1 — 2/i + A, 2 — 2/i + 2A; — , 

g (A,2) = _ m) 2h-A + 0{r l h {535) 
/(A,2) = + 0(r lj ) 

where gi^' 2 ^ and /i^' 2 ^ can also be explicitly written as a sum over two hyper geometric functions. 
As we send £ — >■ oo we recover that 

(#i 2 - (V • Z 1 )(V • Z 2 )) 2 - ^(F • Z X )\V ■ Z 2 f 
G\ >\P X ,P 2 ,Z X ,Z 2 ) iv . Pi)2h -A iv . P2) 2h-A { _ 2Pl .p 2) A . ( 5 - 36 ) 

which is again consistent with the formulas given above. 



6 Numerical results for stress tensors 

The numerical analysis of equation (|5.2|) proceeds largely as for the scalar two-point function, see 
subsection 14. II In particular, we again translate the constraints of crossing symmetry to an infinite 
vector of derivatives at £ = 1 and apply a linear functional in order to exclude certain spectra, using 
the same numerical methods as described above. Notice that the Ward identities (|5.ip can be used 
to express derivatives of / and g in terms of derivatives of h. We therefore do not need to include 
more than the zeroth derivative for the / and g components if we include many derivatives of the h 
component. There is again no guarantee that the coefficients of the conformal blocks are positive 
in the bulk channel. Just as before we will therefore have to assume this condition of positivity in 
order to obtain any bounds. 
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6.1 Bound on the bulk gap 

In order to turn equation (I5.2p into a useful equation to constrain conformal field theories we have 
to decide which parameters we are going to vary. In previous computations of this sort the canonical 
parameter was always the dimension of the external field but for the stress tensor this dimension 
is of course fixed to be d. In our first analysis we instead chose to vary the dimension of the lowest 
spin 2 boundary block which we denote as A( 2 ) ■ We then obtained an upper bound for the lowest 
bulk operator dimension as a function of A( 2 ) which we plotted as the upper curve in figure [TT] 

We may rephrase this result by saying that the upper curve in figure [11] informs us that the 
crossing symmetry equation (|5.2p can only be satisfied if there is at least one "critical" bulk operator 
with a scaling dimension somewhere below the curve. We can however subsequently ask whether 
this "critical" operator really could be sitting anywhere below the curve (and above the unitarity 
bound At, u ik > 1/2). In fact it turns out that the region where such an operator has to appear can 
be constrained even further: we can limit it to the shaded region in figure [TTJ We conclude that for 
every A( 2 ) there has to be at least one bulk operator somewhere within this region. (There could 
in addition be other operators, for example somewhere in the white "band" or multiple operators 
in the shaded region, but none of this modifies the validity of our claim.) 




2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0 



A( 2 ) 

Figure 11: Bounds for the energy momentum tensor two-point function in three spacetime dimen- 
sions. The upper curve is the upper bound A^uik for the first bulk operator as a function of the gap 
A( 2 ) for the first spin 2 boundary operator. The other lines denote further constraints for such a 
bulk operator, to the extend that for every A( 2 ) there has to be at least one bulk scalar somewhere 
in the shaded region. 

In figure [11] we assumed that the vector block Wcis not present in the boundary OPE of T^y. 
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Upon repeating the analysis with a vector block we obtained exactly the same curves for A( 2 ) > 3 
(up to small deviations due to the finite numerical precision), whereas for A( 2 ) < 3 we would not 
be able to bound the bulk gap at all. The latter phenomenon has an easy explanation: the bulk 
identity operator can be decomposed in the boundary channel into the scalar block, the vector block 
and an infinite series of spin 2 blocks starting with A( 2 ) = 3. For A( 2 ) < 3 and with the vector 
block present it is therefore possible to have an infinite gap in the bulk (i.e. no bulk operators apart 
from the identity) and so Abuik cannot be bounded. This is reminiscent of the "trivial" solution 
for the scalar two-point function discussed in appendix IB .51 which we found numerically in section 
1421 

The curves shown in figure have several "bumps" and other features whose origins are unfortu- 
nately unclear to us. For example, we were unable to find specific solutions of crossing symmetry 
that reflect the existence of these bumps. It would be interesting to see if such solutions exist and 
whether a conformal field theory is associated to them. 

6.2 Bound on OPE coefficients in the three-dimensional Ising model 

In subsection 14.3.21 we discussed how to bound OPE coefficients in the conformal block decompo- 
sition. Here we repeat the same procedure for the two-point function of the stress tensor. We will 
again bound the coefficient of the boundary operator T^d which in equation (|5.2p corresponds to 
the coefficient fj3 Q s of the scalar block in the boundary channel. In addition, we decided to focus 
our attention on the three-dimensional Ising model. In particular, we have assumed that the bulk 
spectrum consists of operators with dimensions equal to 1.41, 3.84, and any operator with a scaling 
dimension greater than 4.6. We then obtain an upper bound on /j% ^ as a function of the unknown 
scaling dimension A( 2 ) of the lowest spin two operator in the boundary channel. We assumed that 
no vector operator was present in the boundary channel. Our results are plotted in figure [T2l 

We find a rather surprising plateau for A( 2 ) between approximately 2.9 and 3.2 where ///L ~ 
11.5. From the results in appendix [O we find that /x?qs = 4 in two dimensions and that //? \ = 
640/e + O(e°) in 4 — e dimensions so at the very least our estimate appears to have the right order of 
magnitude. It would be interesting to compute the dimension of the first spin 2 operator appearing 
in the boundary channel in the epsilon expansion, since it is natural to expect that the Ising model 
lies at one of the corners of this plateau. 
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2.5 3.0 3.5 

A (2) 

Figure 12: Bounds for the coefficient of the scalar boundary block in the two-point function of the 
stress tensor as a function of the gap A( 2 ) in the spin 2 boundary dimensions. 

7 Conclusions 

In this paper we have begun to explore the constraining power of crossing symmetry for BCFTs 
in general spacetime dimension. After discussing the basic setup in section [21 we illustrated the 
relative simplicity of the "boundary bootstrap" in section [3] where we found exact solutions with 
at most two blocks in each channel. We have then applied the linear programming methods of [I] 
to the boundary crossing symmetry equations for both scalar operators and stress tensors. With 
our assumption of positivity for the bulk expansion coefficients, we have demonstrated that these 
methods can be useful in the BCFT setup as well and that they lead to interesting universal 
bounds on scaling dimensions and OPE coefficients. Several of our results warrant a more detailed 
theoretical investigation. For example, the bound on the second boundary operator in the special 
transition and the T^d OPE coefficient in the extraordinary transition should be compared with 
computations in the epsilon expansion. Similarly, our numerical results of section I4T21 indicate that 
the bulk-to-boundary OPE always has to be singular, a result that should be put on a more solid 
theoretical footing. Finally, our results for the stress tensor are rather mysterious and certainly 
call for further investigations, beginning with the one-loop anomalous dimension of the spin two 
boundary operator in the extraordinary transition. 

It is unfortunate that the distinct "kinks" of |13] appear not to be generically present in the 
BCFT bounds. We emphasize that this (negative) result is completely independent from our 
positivity assumption, indeed in d = 2 we see no kink but we know that the exact result does 
exhibit positivity. It would be interesting to see if there is another solution to crossing symmetry 
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"standing in the way" and thereby preventing us from obtaining such a kink. More generally our 
results are a reflection of the fact that there are currently no deep understanding of why and when 
such kinks will appear. It would of course be very interesting to understand this phenomenon 
better. We hope hat our numerical results (as well as the analytical results for the minimal models 
of appendix IB .![) will be helpful in further investigations. 

The weakest point of our analysis is admittedly the assumption of positivity for the bulk ex- 
pansion coefficients. While we have presented strong evidence that it is satisfied for the special and 
extraordinary Ising BCFTs, it would be desirable to find a proof. A possible approach would be to 
derive rigorous inequalties for boundary correlators on the lattice. 

It is clear that the possibilities for further numerical exploration are practically unlimited. This 
paper is a first attempt to investigate the boundary bootstrap with a focus on the three-dimensional 
Ising model, but we feel we have just scratched the surface and that there are many interesting open 
questions. To mention a few, there are many other operator dimensions and OPE coefficients to 
be bounded, one may extend our results to supersymmetric theories, and to spacetime dimensions 
greater than four. Furthermore, the relatively simple form of the conformal blocks makes the 
boundary bootstrap especially suitable for investigations involving tensor operators, a research 
direction that is much more involved for the bulk four-point function in a theory with no boundary. 
Finally there is the prospect to broaden the setup and include conformal defects of all possible 
codimensions. 
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A Scalar conformal blocks 

In this section we will use the method of [37] to obtain the scalar conformal blocks as eigenfunctions 
of the conformal Casimir operator. This procedure can be applied with no major changes to two- 
point functions involving tensor operators, and it was used succesfully in section [5] to decompose 
the two-point function of the stress tensor. 

Bulk channel 

The SO(d + 1, 1) generators are, 

Lab = Paq^-P Bwx , (A.l) 
where P A = (P + , P~ , P 1 , . . . P d ). To obtain the conformal blocks we solve the eigenvalue problem 

L 2 (0 1 (P 1 )0 2 (P 2 )) = -C A , (O 1 (P l )O 2 (P 2 )) , (A.2) 

with L 2 = ^(L% + L ( ^ B ){L^ AB + AB ) and C A ,i = A(A - d) + 1(1 + d- 2), where A and I are 
the dimension and spin of the internal operator. Because of Lorentz invariance no operators with 
spin can ever appear in the bulk conformal block decomposition, hence we set I = in equation 
(POL 

Once the asymptotic behavior of /(£) is given, the conformal block is completely fixed. In the 
£ -> limit the bulk OPE dictates [26], 

„ r i(A 1+ A 2 -A) _ (A 3) 

where Ai and A2 are the dimensions of the external operators. Stripping out this factor /(£) = 
£-2(Ai+A 2 -A)g^ an( j pi U ggi n g i n (jA.2p we obtain a standard hyper geometric equation, 

^(1 + 09"(0 + (c+(a + b+ i)Og'(0 + a^(0 = o , (A.4) 

with, 

a = i(A + Ai- A 2 ), 6= i(A- Ai + A 2 ), c = A 
The conformal block for the bulk channel is then, 

/(£) = r^ (Al+A2 - A) 2^i (^(A + Ax - A 2 ), i(A — Ai + A 2 ), A - 
in perfect agreement with |26j . 



+ (A.5) 

- \ + 1; -e) , (A.6) 
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Boundary channel 

In this channel we consider the restricted conformal group. The SO(d, 1) generators are, 

L ^ = P «QpE- P »Qp-a- ( A - 7 ) 

where P a = (P + , P~ , P 1 , . . . P To obtain the conformal blocks we act with the Casimir 

operator on one of the fields and solve the eigenvalue problem, 

L 2 (0i(Pi)0 2 (P 2 )> = -Ca, o (0i(Pi)0 2 (P 2 )) . (A.8) 

where Ca,i = A(A— d+l)+l{l+d— 3) in this case. For this particular two point function only scalar 
blocks are present, so I = again. However, this is no longer true for operators with indices (see 
subsection I5.4j) . The asymptotic behavior for £ — > oo can be obtained from the bulk-to-boundary 
OPE [26], 

/(0~r A . (a.9) 

Stripping out this factor and plugging in (|A.8|) we obtain another hyper geometric equation. The 
boundary block is, 

/(e) = C\F X (A, A - d - + 1, 2A + 2 - d; - i) , (A.10) 
again in perfect agreement with |26j. 

B Solutions to crossing symmetry for scalar operators 

In this section we discuss a few solutions to the crossing symmetry equations for scalar two-point 
functions. We have: 



where A is the conformal dimension of the operator O. The conformal block decomposition is, 



G (0 = 1 + E A fc«fe/buik(A fc ; = e A ^ ^/bd y (A,; (B.2) 

with Afc and /i^ three-point couplings and a& the coefficient of the one-point function of the fc'th 
operator. 
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B.l Two-dimensional Ising model 

In this section we will decompose several correlators for the two-dimensional Ising model. The 
basic fields of the theory, corresponding to the energy and spin operators, will be denoted by e and 
a respectively and have scaling dimensions A e = 1 and A CT = g respectively. As we discussed in 
section there are three different conformally invariant boundary conditions (or boundary states), 
given in equations (I3.18P and (13.191) . The first two are related by the Z2 symmetry of the theory 
and result in the same two-point function of a. The two remaining possible two-point functions for 
the a field are then [25], 



As we shall see below, the + sign corresponds to the extraordinary transition, i.e. the and |e)) 
Cardy boundary states, whereas the — sign corresponds to the ordinary transition which is the \a)) 
Cardy boundary state. 

The full conformal block decomposition can in principle be obtained from Virasoro represen- 
tation theory. We content ourselves here with a simpler analysis where we expand the correlation 
function in the limits £ — > and £ — > 00 and match the coefficients of the expansion to conformal 
blocks. The bulk block decomposition becomes 

= 1 ± |/buik(l; + ^/buik(4; + i ^/buik(8; £) ± ^^/buik(9, + • • • (B.4) 

The bulk spectrum corresponds to the identity 1 and the energy e operators plus scalar Virasoro 
descendants. For example, we may identify the operator of dimension 4 with L_2^-2l and the 
operator of dimension 9 with a level four descendant of e. The absence of an operator of dimension 
5 is in agreement with the fact that e has a null descendant at level two, so L_2^-2£ is actually an 
SO(2, 2) descendant. 

In the boundary channel we find that: 

The constant term in the + case corresponds to a one-point function of a and therefore the Z2 
symmetry is broken by the boundary conditions. We can thus identify it with the extraordinary 
transition. As an additional check one may verify that the bulk block decomposition agrees with 
the decompositions ()3.19j) and (|3.18p . 
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For completeness, we present the conformal block decomposition for the energy two-point func- 
tion. We have [25] . 

Gt = £ + . (B.6) 



so this expression is valid for both boundary conditions. The decomposition in the bulk channel is, 

OO /r) _ o\ -1 

G - = 1 + E W2n;C). (B.7) 



. n- 2 . 

n=l x ' 

For the boundary expansion we obtain, 

r A -G^ = l + f;(^:^" 1 A dy (2n;0. (B.8) 

n=l ^ ' 

From the expressions above we learn that the coefficients of the conformal blocks are positive both 
in the boundary and in the bulk channels. 

B.2 The unitarity minimal models and their analytic continuation 

Let us now generalize the results of the previous subsection to the whole series of the unitarity 
minimal models. Primary operators in the (m, m + 1) model, m > 3, are labeled by integers (r, s), 
with 1 < r < m — 1, 1 < s < m and the identification (r, s) ~ (m — r, m + 1 — s). Denoting the 
(1, 2) operator by a and the (1, 3) operator by e, the relevant OPE and scaling dimensions are 

ax a = l + e, A a = l~ 3 A e = 2-— L^. (B.9) 

2 2(m + 1) m + 1 

We can eliminate m to find 

A £ = ^(4A CT + 1), (B.10) 

and we will work with A CT rather than m as our independent variable from now on. 

We are after the {era) correlator with the Cardy boundary condition labelled by the identity. 
(Recall that in the Ising model this Cardy state is associated to the extraordinary transition, see 
equation (I3.18P ). This correlator can be obtained as a special case of a result obtained in the 
context of Liouville theory with ZZ boundary conditions [38], where the two-point function 

(V- b/2 (x)V„(y)) (B.ll) 

was evaluated. Here V a (x) denotes the usual Liouville vertex operator with scaling dimension 
A a = a(Q — a) with Q = b + b^ 1 . We will be interested in the case a = —6/2 and b set to the 
minimal model values given by 

C= 1 + 6 Q 2 = 1- 6 (B.12) 
m(m + 1) 
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One may verify that solving this equation for b results in a scaling dimension of V_b/2 which is 
precisely A a given in (|B.9p . 

The two-point function from [38j takes the form: 



G aa (0 = 2 sin (~(1 + 4A CT )) ^+^(1 + ^ +1)/3 
1-2A„ 2 + 2A„ 2-4A„ 1 



X2F1 



(B.13) 



3 ' 3 3 f + 

The boundary conformal block decomposition of this correlation functions contains operators with 
even dimensions, 



,-A. r f M_o. rn 1 , A x\ I, , A„(l + A g ) , 9 ,s , A CT (l + A a ) 2 (2 + 5A g ) 
e G CTCT (0-2sm(-(l + 4A (r )j ^1+ 2(5 _ ^ / bdy (2, £) + 4Q(n _ 4 ^ )(g _ ^ /bd y (4, 



A g (l + A ff ) 2 (20 + 106A. + 35A 2 + 21A 3 ) v 
+ 1008(17 -4A g )(ll-4A g )(5-4A g ) + l^j 

in agreement with the fact that the only boundary block is the identity Virasoro block. In the bulk 
channel we find the identity and the e Virasoro blocks, leading to a decomposition into SO(2, 1) 
blocks with operators of dimension of 1 + 4n and A £ + An with n a non-negative integer. For the 
first few coefficients we find 

A g (l + A g ) A g (l + A g ) 2 (2 + 5A g ) 

GW(0 - 1 + 2(5 _ 4A<t) W4,0 + 40(11 -4A CT )(5-4A (T ) /bulk(8 ' e) + " " " (R15) 

r(l=j^)T(m^)( (l + A g )(2 + 5A g )(-l + 8A g ) , A , 4 >x, 

Jbulk(A e , H /bulk(A £ + 4, 4 j + . 



r(-2A g )r (^±|^) V 6(7 + 4A (T )(5 + 8A (J ) 

Up to the normalization factor 2sin(^(l + 4A g )), the coefficients of the first series are the same as 
those of the boundary identity Virasoro block. Indeed, in either channel these blocks correspond 
to Virasoro descendants of an identity operator. Notice also that the coefficient of the block with 
dimension A e + 4 has a zero precisely when A a = |, reflecting the aforementioned fact that L_2£ 
is actually an 50(2, 2) descendant in the two-dimensional Ising model. (Indeed in the Ising model 
the (1,3) primary is identified with the (2, 1) primary which has a level-two null descendant.) 

Remarkably, the coefficients of the boundary conformal blocks turn out to be positive for < 
Ac- < |@ This implies that we have found a solution of the crossing symmetry equation that is 
consistent with the unitarity requirements for any value of A CT in this interval, given simply by the 
analytic continuation of (|B.13P away from the minimal model values for A a . Of course this does 



9 We have verified this statement to high order and believe that it is generally true although we currently cannot 
offer a rigorous proof. 
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not imply that this correlator can always be embedded in a full-fledged unitary CFT - in fact we 
already know that this is only possible if A CT has one of the minimal model values. 

As we pointed out repeatedly in this paper, unitarity does not require the coefficients of the 
bulk channel conformal blocks to be positive. These coefficients however do turn out to be positive 



symmetry equation ()2.2ip . with the dimension of the first bulk scalar primary e in the a a OPE 
given by (jB.lOh . In the smaller range | < A CT < 1 the bulk expansion satisfies positivity. 

An aside: the four-point function 

The result just found compels us to briefly consider the analogous analytic continuation for the 
mimimal- model four-point function (without boundary). The numerical bounds in that case [5] 
appear to converge to a shape with a "kink" at the Ising model, so at A CT = |, which is followed 
by a straight line that is approximately given bv lBTTOl for \ < A CT < 0.4. (For larger values of A CT 
the numerical analysis becomes less precise.) 

Now, as we show momentarily, one can easily repeat the boundary analysis and construct a 
four-point function {aaaa) with A £ given by IB. 101 This solution appears to precisely saturate the 
numerical bounds of [5j in the range | < A CT < ^, which explains why they cannot be lowered 
by e.g. improving the numerical accuracy, even in between the discrete minimal model points. 
Note of course that there is no complete unitary CFT when one does not precisely sit at the the 
minimal model points, only a solution for this particular four-point function. On the other hand, we 
currently cannot explain why this solution is "extremal" in the sense that it saturates the bounds. 

We can construct this "interpolating" solution by noticing that the a field in the minimal 
models has a null descendant at level two and its correlation functions therefore satisfy the following 
differential equation: 



for the smaller range | < A CT < 1. The lower and upper endpoint of this range are determined by 
the zeroes of the blocks of dimension A e + 4 and A e , respectively. 

In summary, for the range < A CT < | we have found an exact solution to the boundary crossing 




(B.16) 



with 



(I 



Z 1 




(B.17) 
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Our natural conjecture is that the same differential equation is satisfied by a putative four-point 
function that saturates the bound and interpolates between the minimal models. 

When we apply this differential operator to the four-point function of four a fields, which we 
write as usual as 

{(j{x 1 )a(x 2 )a{x 3 )a(x 4 )} = ^—^—^ G{z, z) , (B.18) 

we find a simple hyper geometric equation that we can easily solve. Combining the holomorphic 
and antiholomorphic part the two solutions become: 

G(z, z) = G 1 [z)G 1 (z) + N(A)G 2 (z)G 2 (z) , 
Gi(z) = (1 - z)-\ Fl (1^* -2A, s) = 1 + 0(, 2 ) , (B . 19) 

G 2 (z) = (1 - z)^z^ 2Fl (^±^, 1 + 2A, 4I±*> , ^ = Z (4A +1 )/ 3(1 + 0{z)) . 
Crossing symmetry fixes the relative normalization to be: 

N(A) = g^Z!Z£ (I - r(l + 2A) 2 (- cos (U- + 4irA)) + sin (| (vr + 16ttA))) 

vrr(| + f) 2 

Notice that iV(A) > for < A < 1 but JV(0) = N(l) = 0. 

The conformal block decomposition can be found most easily by first decomposing the holo- 
morphic functions separately in terms of the holomorphic building blocks 

f(P,z) = zP 2 Fi{P,P,2p,z). (B.21) 

For the identity block we then find that 

oo 

Gi(z) = X>(n)/(2n,z), (B.22) 

n=0 

with the first few coefficients given by: 

ci(0) = l, 

A(l + A) 



ci(l) 
d(2) 
ci(3) 



2(5 - 4A) ' 

A(l + A) 2 (l + ^) (B-23) 
20(-ll +4A)(-5 + 4A) ' 
A(l + A) 2 (5 + |A (53 + |A(5 + 3A))) 



252(-17 + 4A)(-11 + 4A)(-5 + 4A) 
and we checked that the first six coefficients are all positive functions for < A < 1. 
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The epsilon block can be decomposed as: 



G 2 (z) = V c 2 (n)f ( ^-tl + 2n, ^ , (B.24) 

n=0 V 6 ' 



with 



c 2 (0) = l, 

(1 + A) + f ) (-1 + 8A) 
C2Uj 3(7 + 4A)(5 + 8A) ' ( R25 ) 

(1 + A) 2 (-35 + ±A (419 + ±A(6315 + 64A(97 + 25A)))) 
~ 18(7 + 4A)(13 + 4A)(11 + 8A)(17 + 8A) ' 

We again checked by hand that the first eight coefficients are positive for | < A < 1. They however 
also all have simple zeroes for some A < g, starting with 02(1) at A = |. This first zero corresponds 
to the decoupling of a spin 2, dimension 3 operator (as well as an infinite number of its Virasoro 
descendants) which we can again identify with the level two null descendant of the epsilon operator 
in the two-dimensional Ising model. The fact that our putative solution extends to A = 1 implies 
that it extrapolates beyond the accumulation points of the minimal models at A = ^. 

Notice that for A CT > ^ it is natural to conjecture that the bound will be saturated instead 
by the free boson where (in the notation of this section) A t = 4A CT . On the other hand, for 
A CT < | our putative solution becomes invalid precisely because the coefficient of the conformal 
block corresponding to L_2£ becomes negative. It would be very interesting to find a solution 
saturating the numerical bound for A a < g which would would allow one to investigate the precise 
transition at the "corner" corresponding to the Ising model point and the role of L_2£ in this 
transition. Interestingly, the (approximate) decoupling of the first irrelevant spin 2 operator was 
observed numerically in three dimensions as well [13]. 

B.3 (</> 2 </> 2 ) correlator 

In this section we will decompose (4> 2 4> 2 ) i n f ree field theory. This expansion complements the order 
e expression for the scalar two-point function of section [3j The <p 2 two-point function is, 



C^lifAV") +^" 2 , (B.26) 



A + J 2 

where the plus/minus sign corresponds to Neumann/Dirichlet boundary conditions, and N is the 

number of scalars. The conformal block expansion in the bulk channel is 

oo 

= 1 + \ a< p2f hnlk (d -2-0 + Yl Aa ^,n/buik(2d - 4 + 2n; , (B.27) 



n=0 
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with 

Aa,02 = ±2 , 

_ ((-l) n 2 c! r(^i)r(i(i + n - 1) + 4iV^r(d + n - 2)) r(d + n - 2)r(|d + n - 4) 
° 04 ' n ~ 80Fr((i - 2)2r(n + ljf(p + 2n - 4) ' 

(B.28) 

The Neumann expansion exhibits positivity while the Dirichlet case has one negative coefficient. 
In the boundary channel we have, 

AT 00 

C d+2 G%^ 2 = y + E ~ 2 + 2n > > (B.29) 

n=0 

with 

u*-h+fi r(^ + n)r(id + n-l)r(d + 2n-3) 

^ n -(l±^0j (2n)! r( l d _ 1)r(d _ 2)r( *t^=3) ' [ ^ 6U) 

with positivity in both cases. 

B.4 The extraordinary transition 

There is no extraordinary transition in 4 dimensions since the conformally invariant one-point 
function of a free field is not compatible with its equation of motion. In 4 — e dimensions the 
equation of motion is however modified to: 

U4> = ^ 3 (B.31) 

with A* = 48n 2 e/(N + 8) and N = 1 in our case. On the half-space this equation admits the 
solution: 

12 1 

which to leading order is consistent with boundary conformal invariance. This solution is our 
starting point for the analysis of the extraordinary transition in the e expansion. 

Let us compute the two-point function of the scalar field (p. We may shift the field cp by its 
classical one-point function, 

<Kx) = (<f>(x))+x(x) (B.33) 

and find the propagator G(x,y) = (x{ x )x{y)) by solving the linearized equation of motion around 
this solution, 

(n - ^) G(x, y) = 5 d (x - y) (B.34) 



(<Kx)) = \l--i (B-32) 
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The solution compatible with the boundary conditions at x d = takes the form: 



G(x,y) 



{2x d )(2y d 



h^-f 1 (i G0(o ) 



I ( £ \ (B - 35) 

with £ = (x — y) 2 /(Ax d y d ), as before. On the first line we recognize the familiar form of a scalar two- 
point function for a CFT with a boundary. Taking the limit £ — > in ()B.35|) we see that the properly 
normalized operator is actually 2irx rather than and similarly 27T(fi rather than (j). We will 
henceforth work with these rescaled operators. This implies that from now on (4>(x)} = 3/(^x d ) 
and we can drop the Att 2 on the first line of (jB.35 j) . 

We will now expand the two-point function of (j> in conformal blocks. It is important to note 
that OPE statements always refer to full correlation functions, i.e. including any disconnected 
contributions. In our case the disconnected part {4>{x)) {<p(y)) is of order 1/e which makes it the 
leading-order term. Our first task is thus to decompose (<p(x)) {(p(y)) in conformal blocks. In the 
boundary channel we of course find precisely the block corresponding to the identity operator and 
nothing else. In the bulk channel we find: 

<«*»<«»» = pl|7) = (W) rl ( A " lk(2 ' + |w A - (2 + 2 "' {) ) (R36) 

Interestingly, the product of two one-point functions decomposes into an infinite set of bulk blocks 
with dimensions given by the even integers and with positive coefficients. However, as expected for 
a totally disconnected correlator, the bulk identity operator is missing at this order. 

At the next order we should take into account that the one-point function of (ft a priori has 
subleading corrections, 

(</>(x)) = (l + e a + | log(2x d )) (B.37) 

with an unknown coefficient a and with the logarithm originating from the correction to the scaling 
dimension of <p in 4 — e dimensions. The full two-point function to order e° becomes: 

<«»)*)> - ( ^4 V)A . c^"«> (b38) 

G ext (0 = -(1 + 2ea)£ - 18£log(£) + G°(£) 
e 

with A (A = 1 — e/2 the free-field dimension of <j> in 4 — e dimensions. 

In the boundary channel the conformal block decomposition of this corrected correlator is again 
straightforward. The corrections to the disconnected part of course simply become corrections to 
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the boundary identity block, whilst for the connected part we find that: 

r 1 G°(0 = ^/bd y (4;0 (B.39) 

so we find a single boundary block of dimension d = 4. This is completely as expected. In particular, 
the existence of a gap of size d was an essential assumption in the numerical bootstrap for the bulk 
bounds. 

In the bulk channel we find subleading corrections to the infinite series of blocks in (IB.36j) but 
no new blocks. The first few terms take the form: 

G ext (0 = 1 + (™ + 11 + 72a) / bulk (2 - ~e; £ 

+ - 12 + 72a) /bulk (4,0 + " 18 + 2 ! ■A-n,(<> + 2r.O (B. 10) 



+ (i + 2i) ( - 241 + 144a) ) /bulk ( 8+ T^) + 



where it is understood that the blocks are evaluated in 4 — e spacetime dimensions. The order 
1/e terms in (1B.40P of course coincide with (IB. 36ft . The identity operator is now present in the 
bulk channel, and the dimension of the next operator (which is 2 — |e) is precisely the one-loop 
dimension of (f) 2 in the epsilon expansion. It would be interesting to compute a so we can get an 
idea of positivity of the coefficients for e = 1. 

B.5 A trivial solution 

A particularly simple solution of (|2.2ip is obtained by assuming that the bulk channel only contains 
the identity operator, so all the non-trivial one-point functions are set to zero. In that case there 
is effectively no boundary at all and the two-point function is just (x\ — X2)~ 2A . This two-point 
function still has a boundary conformal block decomposition of the form: 

00 

r A = J>mAdy(A + m;0 (B.41) 



with 



m=0 



2^(A) m (A - I + l) m /2(A - V" + m)_ m / 2 m even 

1 (A) m (A-f + l)( m _i)/2(A-^fi + m) (1 _ m) /2 m odd 



(B.42) 



2 m m! ' 

All the coefficients are positive for A greater than the unitarity bound and the boundary spectrum 
begins with an operator of dimension A. 
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B.6 Generalized free field theory 

As a simple generalization of the free field theory result we define generalized free field (or gff) 

two-point functions in the presence of a boundary as: 

1 1 



<0(xi)0(x 2 )) - , 



(xi-a;2) 2A ((xi -x 2 ) 2 +4xf^) 

r A c&(0 g±(o = i±' 



(B.43) 



The conformal block decomposition in the bulk takes the form 

» (-l) n (A) n (2A-|+2n) 

n=0 V 2 / — n 

which has the expected 'double trace' infinite operator spectrum and coefficients with alternating 
signs. On the boundary we find that: 

~ (AW (A - M + 2n) 

2 4 

r A,- m ^ (A) 2n+1 (A-^ + 2n)_ n r _ ^ 

* = g 2 2 »(2n + D! (A - I + n + l)_ w ^ + + ^ 

and we find two 'single trace' operator spectra on the boundary, both with positive coefficients. 

B.7 O(N) model at large N 

For the the O(N) model with Neumann boundary conditions the scalar two-point function is given 
by [26], 

1 x^- 1 



G 0(N) = [j-^) (1 + 20- (B.46) 
The bulk channel expansion is, 

oo 

G (N) =1 + J2 Aan/buik(2n + 2; £) , (B.47) 

n=0 

with 

A2n (rf 2 - 4rf(n + 2) + 8(1 + n) 2 + 4)r(l -\d + n)T(2 -\d + nf 

Aa n = ( — 1) — 2 1 • (B.48) 

4r(2 - \d) r(n + 2)r(2 - id + 2n) 

As in all the expansions with Neumann boundary conditions studied in this appendix, the bulk 

channel coefficients are positive. Finally, the boundary channel expansion is, 

C ¥+l G 0{N ) = 2/bd y (d - 3; . (B.49) 
It is somewhat unexpected that in this channel we have a single block. 
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C Conformal block decompositions for T^ v 

In this appendix we present a few explicit examples of conformal block decompositions of the form 
(|5.2p for the two-point function of the stress tensor. 



C.l Two bulk dimensions 

The conformal block decomposition of the stress tensor two-point function in two dimensions is a bit 
subtle, see for details. First of all, the residual Virasoro symmetry plus the absence of energy 
flow across the boundary completely determines the two-point function. Furthermore, the number 
of independent tensor structures decreases to two and the two functions /(£) and g{£) have to be 
replaced with the single function 2£/(£) + (1 + With our unit normalization we find that 

the resulting two-point function is given precisely by the boundary scalar block, with a coefficient 
that is equal to 4. In the bulk channel we find the identity plus a single block of dimension 4 with 
unit coefficient. 

C.2 Free field theory for general d 

The two-point function of the stress tensor in free field theory for d > 2 decomposes into infinitely 
many blocks in either channel. Without presenting all the formulas, we have presented the first few 
operators and their associated coefficients in both the bulk and the boundary channel in the tables. 
Notice that the coefficients in the bulk channel are not positive for either boundary condition. 

C.3 Extraordinary transition 

In this subsection we compute the two-point function of the stress tensor in the extraordinary 
transition. 

The classical stress tensor for the A^ 4 theory with a curvature coupling z takes the form: 

T^(x) = ff 2 z - d^d^ + 2z (/tipW + g^ + Q - 2^ d p m - 2z<j>U^j ) 

(C.l) 

One may easily verify that it is traceless in d = 4 for z = 1/12 which therefore corresponds to the 
conformally coupled scalar. We will henceforth use z = 1/12. In that case T^ v is unit normalized 

H 2 i 

in free field theory, more precisely (T^(x)T p(7 (y)} = { ^ )i{ Jy {yd)i provided (0(a#(y)) = ^=^r- 

The correlation functions of the scalar <j) were computed to leading order in subsection IB.4I 
Upon substituting the solution (4>(x)) = 2>/{^fex d ) in (IC.ip we find that the one-point function of 
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A 







1 


d-2 


_L / (-2+^(1+^ 

* ^ 4(-i+d) ; 


2d 


+i 


2d + 2 
2d + 4 
2d + 6 


f(-2+d)<A 

V (2+3rf) ^ 
/(-2+d)d 2 (l+d)\ 
+ ^ 6(2+d)(4+3d) ^ 
/(-2+d) ( i 2 (l+d)(2+d)\ 
^ 18(8+3d)(10+3(i) J 


2d + 2m 





Table 2: Bulk conformal block decomposition of the two-point function of the stress tensor in 
free field theory. The first block corresponds to the identity operator and its coefficient sets the 
overall normalization. The plus/minus sign corresponds to the special/ordinary transition, i.e. 
Neumann/Dirichlet boundary conditions. 



A 


I 




d 





2d 

(-1+d) 


d 


2 


2 1 - M (1±1) 


d + 2 
d + A 


2 
2 


2- 2 - 2d d(-l+d)(2+d) 
2- & - 2d d{-l+d)2+d) 2 (4+d) 


3(7+d) 


d + 2m 


2 





Table 3: Boundary conformal block decomposition of the two-point function of the stress tensor in 
free field theory. 

T^ v vanishes, in agreement with the requirements of boundary conformal invariance. At the next 
order we substitute (ft(x) = (</>(x)) + xi x ) an d expand in e to find an expression of the form: 

T pu (x) = ^-T^ d ,d x ] X (x) + ... (C.2) 

where T t _ l u[x d ,d x ] is a linear differential operator which explicitly depends on x d . To leading order 
we therefore obtain that 

(T^(x)T pa (y)) = -T^[x d ,d x ]T p(T {y d ,d y ](x(x)x(y)) (C.3) 

We can now substitute the solution G(x,y) = {x{ x )x{y))i which is equation (|B.35P without the 
factor of 4-7T 2 , work out the action of the differential operators T and collect various terms to 
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eventually find a two-point function of the form: 



(T^(x)T p(7 (y)) 



(C.4) 



(4£)4(z^(^)4 

where the tensor structures H\2 and Qu are defined (in the projective cone notation) in (|5.7|) and 
(|5T8j) and 

64£ 2 (2 + 50 ,_ ext/ ^ 64£ 3 (1 + 5£ + 10£^ 



/ ext (0 



16f 



5 ext (0 



/i ext (0 



(C.5) 



e(l + 3 e(l + 4 e(l + 5 

Upon comparing (lC.4j) with the last equation in (I5.17D we see that this correlation function has ex- 
actly the right tensor structure to be consistent with boundary conformal invariance. Furthermore, 
the functions (/ ext , g ext , h ext ) also satisfy the Ward identities (j5. 19|) . These are rather non-trivial 
checks of our result. 

The conformal block decomposition of ()C4|) turns out to be remarkably simple. In the boundary 
we find only a scalar block (which must have dimension d by the Ward identities) with coefficient 
640/e. In the bulk we find three blocks, 

// ext (0\ 



,ext (e) 



\/i ext (0/ 



160 

e 



//buik(2;0\ 

Sbuik(2;£) 
Wuik(2;£)/ 



+ 



480 



//buik(4;£)\ 



flt>uik(4;f) 
V/ibuik(4;e)/ 



+ 



320 



/7buik(6;0\ 



5buik(6;C) 
V^buik(6;0/ 



(C.6) 



all with positive coefficients. Notice that the identity operator is absent at this order. 

Closer inspection of (|C.6j) leads to a subtlety that we would like to clarify. We easily identify 
the bulk block with dimension 2 as the operator cj) 2 . It appears in the TT OPE with an order one 
coefficient and its one-point function is {4> 2 ) = ((f)) 2 ~ e _1 so altogether it appears at the right 
order in e. The counting for the operator of dimension 4 is however a bit different. The only scalar 
primary of that dimension is <f> but its one-point function is of order e~ 2 . Our result can therefore 
only be consistent if cp 4 appears in the stress tensor OPE only at order e. It is in fact easy to see 
that the leading-order Feynman diagram for the (TT0 4 ) tree-point function (which would be of 
order e°) has to vanish. This is because it factorizes into a product of two Feynman diagrams that 
each correspond to the {T(j) 2 ) two-point function, which in turn vanishes by conformal invariance. 
This is also consistent with the fact that no dimension 4 block appears in the bulk conformal block 
decomposition of the stress tensor two-point function in free-field theory, cf. table [2j From these 
tables we may also deduce that a similar cancellation should occur for the dimension 6 operator. 
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